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ABSTRACT 

Detailed chemical abundances for twenty one elements are presented for four red giants in the 
anomalous outer halo globular cluster Palomar 1 (i?GC = 17-2 kpc, Z = 3.6 kpc) using high-resolution 
{R = 36000) spectra from the High Dispersion Spectrograph (HDS) on the Subaru Telescope. Pal 1 has 
long been considered unusual because of its low surface brightness, sparse red giant branch, young age, 
and its possible association with two extragalactic streams of stars — this paper shows that its chemistry 
further confirms its unusual nature. The mean metallicity of the four stars, [Fe/H] — — 0.60± 0.01, is 
high for a globular cluster so far from the Galactic center, but is low for a typical open cluster. The 
[a/Fe] ratios, though in agreement with the Galactic stars within the la errors, agree best with the 
lower values in dwarf galaxies. No signs of the Na/0 anticorrelation are detected in Pal 1, though Na 
appears to be marginally high in all four stars. Pal I's neutron capture elements are also unusual: its 
high [Ba/Y] ratio agrees best with dwarf galaxies, implying an excess of second-peak over first-peak 
s-process elements, while its [Eu/a] and [Ba/Eu] ratios show that Pal I's contributions from the r- 
process must have differed in some way from normal Galactic stars. Therefore, Pal 1 is chemically 
unusual, as well in its other properties. Pal 1 shares some of its unusual abundance characteristics 
with the young clusters associated with the Sagittarius dwarf galaxy remnant and the intermediate- 
age LMC clusters, and could be chemically associated with the Canis Majoris overdensity; however it 
does not seem to be similar to the Monoceros/Galactic Anticenter Stellar Stream. 
Subject headings: galaxies: dwarf — globular clusters: general — globular clusters: individual(Pal 1) 



1. INTRODUCTION 



Ever since its discovery by lAbelll (|1955| ) , Palomar 1 
(Pal 1) has been tentatively classified as a globular clus- 
ter (GC), primarily because of its location high above 
the Galactic plane {Z = 3.6 kpc; Harris 1996, 2010 
edition). However, Pal 1 has several anomalous char- 
acteristics that have cast doubt upon this classifica- 
tion. For example. Pal 1 appears to be unusually young 
for an outer halo GC. Using standard isochrone fits to 
Pal I's ground-ba.sed co lor-magnitude diagram (CMD), 
iRosenberg et al.l (11998a) es timated an age between 6.3 
and 8 Gyr, whilelSarajcdin i et al. | (J2007) found an age 
between 4 and 6 Gyr with HST photometry. This F606W 
(^V ) F814W (~I) HST CM D of the central field of Pal 1 
from lSaraiedini et al.ll2007l (see Figure [Ij shows the main 
sequence turn-off around F606W w 19.4 {V w 19.6), 
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and a possible detection of a red horizontal branch near 
F606W « 16.4 {V « 16.6). Thus, Ay^B^o « 3, 
making the cl uster younger than a bout 7 Gyr (with 
the method of IChabover et al] 119961) . These estimates 
make Pal 1 one of the youngest GCs, along with Ter 
7, Pal 12, and Whiting 1, whose age s range from 7 — 9 
Gyr (e.g. .B uonanno ct al. 1998; Sa laris fc Weiss! 120021: 
ICarraro et al.ll2007l : IDotter et al.ll2010[ ). 

Pal 1 also lacks luminous, evolved stars. The CMD in 
Figure [T] shows a sparsely populated red giant branch 
and a barely detectable horizontal branch. Further- 
more, near-infrared medium-resolution spectroscopy of 
the Ca II triplet (CaT) lines in four of these red giants has 
shown that Pal 1 is fairly metal-rich for an outer halo GC, 
with an average [Fe/H] = —0.6 (on the Zinn-West scale, 
using the standa rd Galactic CaT line str ength to metal- 
ficitv conversion: [Rosenberg et al.|[l998bD . The main se- 
quence is also poorly populated, and therefore the entire 
cluster has an extremely lo w total luminosity for a GC 
(My = -2.52: lHam3[T99l 2010 edition). 

Clearly Pal 1 does not conform to the characteristics 
of traditional Galactic GCs, suggesting either that Pal 1 
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Fig. 1. — An HST color-magnitude diagram (roughly V, V-I) 
for the central field of Pal 1 from Sarajcdini et al. (2007). The 
observed targets are circled. Also shown are fitted isochrones from 
the Dartmouth Stellar Evolution Database (Dottcr et al. 2008) for 
ages of 4, 5, and 6 Gyr. Values of [Fe/H] = -0.6, [o/Fe] = 0.0, 
(m-M)v = 16.27, and E{B-V) = 0.20 are adopted, as discussed 
in the text 

is not a GC or that it is not Galactic. If Pal 1 has been 
misclassified as a GC, another possibility is that it could 
be an open cluster (OC) — this option seems unlikely due 
to the combination of its position above the Galactic 
plane , its high concentration parameter (c = 2.57: lHarrisl 
119961 2010 edition), and its age (although Pal I's age 
is comparable to that of th e oldest known OCs, e.g. 
iVandenBerg fc StetsonI 120041 ) . However, Pal I's bright- 
ness and age place it in greater ag reement with the OC s 
than with the GCs (see Figure 3 in lCarretta et al.ll2010[ ). 
The ambiguity surrounding the OC vs. GC classifica- 
tion leads to a need for a more rigorous distinction be- 
tween the two. From o bservations of many Galactic GCs, 
iCarretta et al.l ()2010f ) defined a GC as a cluster that 
shows signs of a second generation of stars through a 
Na/0 anticorrelation. Under this definition, if Pal 1 does 
not show signs of a Na/0 anticorrelation, then it cannot 
be considered a classical GCo 

Alternatively, Pal 1 could have originated in a dwarf 
galaxy and been accreted by the Milky Way during 
a merger. The latter situation is plausible, as sev- 
eral other Galactic GCs, notably the young clusters 
listed above, have been linked to dwarf galaxies. M54, 
Arp 2, Ter 7, and Ter 8 lie on a stream associated 
with th e Sagittarius dwarf spheroidal (Sgr dSph) galaxy 
(|Da Co sta & Armandrofi 1995); Pal 12 has been kine- 
matically lin ked to the Sgr dSp h through proper mo- 
tion studies (jDinescu et al.l[2000l ): Whiting 1 is under- 
going tidal stripping and app ears to lie in a stre am of 
M giants from the Sgr dSph (jCarraro et al.ll2007| ): and 

^ In the past, Pal I's unusual characteristics have led to a 
designation as a transitional cluster, i.e. an object that resem- 
bles something bet ween classical globular and open clusters, e.g., 
lOrtolani eTalT frogs' ) . 



Rup 106 has been tentatively associated with the Magel- 
lanic Cl ouds because of its l ocation along the Magellanic 
stream ()Lin &: Richeii 119921 ) . Pal 1 does appear to lie 
along several stellar streams that may be extragalactic 
but have not yet been associated with any galaxies. One 
such stream, known both as the Monoceros stream and as 
the Galactic Anticenter Stellar Stream (G ASS), appears 
to be from a disrupted dwarf satellite (e.g. iSollima et alj 
120111 ): Pal I's position, velocity, an d metallicity mak e 
it a possible member of this stream (jCrane et al.l 12003^. 
The Canis Major overdensity. which seems to be distinct 
from the GASS (|Chou et al.ll2010bD . has also been as- 
sociated with Pal 1 , again because of the cluster's co- 
incidental location ([Martin et al.ll2004D . Finally, Pal 1 
has also been tentative ly linked to the "Orphan Stream" 
(jBelokurov et al.ll2007| ). though its radi al velocity and 
metallicity make it an unlikely member (jNewberg et al.l 
l2010f l. Though no obvious association with a particular 
stream has been confirmed, tidal tails around Pal 1 have 
recently been discovered, suggesting that the cluster it- 
self is being dissipated in the tidal field of the Galaxy 
(jNiederste-Ostholt et alll2010( ). 

In the absence of an obvious host galaxy, detailed 
chemical abundances can be used to examine if Pal 1 has 
an extragalac t ic ori g in. It has been wel l established (e.g. 
iTolstov et all 120091 iVenn et al.l I2004D that metal-rich 
stars in Local Group dwarf galaxies have distinctly dif- 
ferent chemical abundances from stars in the Milky Way, 
at a given metallicity (the differences are not so obvi- 
ous by [Fe/H] ^ —2). These chemical variations between 
the environments are due to different star formation his- 
tories and efficiencies, which are linked to the mass of 
the galaxy and the properties of its environment. There- 
fore, the differences are not just between dwarf galaxies 
and normal galaxies — in fact, dwarf galaxies have unique 
chemical evolutions, meaning that it is pos sible to chem- 
ically link a star to its host galaxy (e.g . ITolstov et al.l 
120091 iFreeman fc Bland- HawthornI 120021) . These chemi- 
cal differences are not jus t seen in field stars, they are also 
seen in GCs; for example. iPritzl et al. (2005!) showed that 
most Galactic GCs have similar chemical abundances to 
Galactic field stars at the same metallicity, with Pal 12, 
Ter 7, and Rup 106 as notable exceptions. 

Detailed, high-resolution chemical abundances have 
been used to trace the origin of sev eral g lobular clus- 
ters. In particular Pal 12 (Cohen '2004) and Ter 7 
(jTautvaisiene et al.l 12004 |Sbordonc ct al. 20051) have 
[a/Fe] ratios that follow the Sgr dSph field star abun- 
dance trends and are clearly separated from Galactic field 
stars, supporting the hypothesis that they were captured 
from the Sgr dSph. The more metal-poor Sgr clusters 
Arp 2 and Ter 8 . however, are not chemically distinct 
from the Galaxy (jMottini et al.|[2008[) , which shows that 
the [a/Fe] deviation can be metallicity- (and presumably 
age-) dependent. If Pal 1 has an extragalactic origin, 
then its metallicity of [Fe/H] « —0.6 puts it in a regime 
where [a/Fe] and other chemical abundance ratios should 
be useful in determining its origin. 

Here we present a detailed analysis of high-resolution 
spectra for four re d giant stars in Pal 1 , taken with the 
Subaru Telescope. iMonaco et al.l (|2010[ ) have completed 
a chemical analysis of one star in Pal 1 (Pal l-I) , also with 
spectra taken at Subaru, but with independent observa- 
tions; they find that their single star shows similar abun- 
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dance patterns to Galactic open clusters. iSaviane et alj 
(|201Clf) further show that its chemistry implies that Pal 
1 could be associated with the Canis Major overdensity. 
With four stars, our analysis allows us to determine the 
chemical abundance ratios with higher precision, and in- 
vestigate any star-to-star variations in the abundances 
that could be linked with chemical or cluster evolution 
effects. Section [2] outlines the observations and the data 
reduction, while Section[3]examines our method for mea- 
suring equivalent widths. In Section |4] we discuss the 
model atmospheres, particularly how the atmospheric 
parameters are derived. Section [5] presents the methods 
for element abundance determinations, while Section |B] 
examines the significance of these results with respect to 
Galactic stars and clusters. Finally, Section [7] examines 
the origins of Pal 1. 

2. OBSERVATIONS AND DATA REDUCTION 

Four of our targe t stars are the same stars observed 
by iRosenberg et al.l (|1998b( ): bright, red giant branch 
(RGB) stars, with radial velocities that make them likely 
cluster members. We refer to these stars as Pal l-I, -II, 
-III, and -IV, following the naming convention of Rosen- 
berg et al. The bright object located at the center of the 
cluster. Pal 1-C, was also included as a target, although 
the membership of this object is not as certain as the 
others. The locations of the se stars on a CMP an d in 
an image of the cluster ffrom lSaraiedini et al.ll2Q07l ) are 
shown in Figures [1] and [2] In addition to the Pal 1 stars 
we include a n analysis of a "st andard" star, M67-141 (as 
identified bv iFagerholml llQOa ) . This star resides in the 
metal-rich open cluster M67, and has previo us detailed, 
high r esol ution chemical abundances from lYong et aD 
(|2005D and lPandnoet al.l (|2010[ ). Table [H shows the po- 
sitions and magnitudes of the target stars. The Pal 1 
positions and V and I magnitudes (in the Johnson sys- 
tem) are from Sarajedini et al. (2002), while M67-141's 
position and V a nd I magn i tudes ( in the K r on-Co usins 
system) are fro r nlHog et all ()2000t ). iSandersI ()1977f ). and 
iJanes fc SmithI ()1984[) . respectivelv. All of the K magni- 
tudes are from the Two Micron All-Sky Survey (2MASS) 
Point Source Catalogo 

These stars were observed on several runs in 2006 Jan- 
uary and 2007 January using the H igh Dispersion Spec- 
trograph (HDS jNoguchi et al.ll2002D on the 8.2 m Subaru 
Telescope. The default grating was centered at 5500 A 
with a slit length of 5.6 arcsec for all stars and a slit 
width of 1 arcsec for the Pal 1 stars and 0.6 arcsec for 
M67-141, leading to spectral resolutions of i? = 36000 
and R — 60000, respectively. Table [2] shows a list of 
the dates, total exposure times, and signal-to-noise ratios 
(SNR) for these observations. We note that the position 
of the Moon was not preferable for observations of the 
Pal 1 stars, as the primary ta rget of this obser ving run 
was the Sextans dwarf galaxy (|Aoki et al.ll2009D . 

The data were reduced using the basic Subaru HDS 
pipeline, with standard IRAFq routines to remove the 
bias and fiat field, and with additional cosmic ray re- 
moval. The stars were fairly low in the sky during the 

^ pittp : / /tdc-vnm . harvard . edu/catalogs/tmpsc . html | 
"* IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 




Fig. 2. — An HST F814W image of Pal 1, with our targets shown 
in squares. The field of view is approximately 1.'7 X 1.'7. North 
is up, and west is to the left. Note that we do not perform a full 
analysis of Pal 1-IV due to poor sky conditions. 



observations and the seeing conditions were poor, mean- 
ing the stars observed in 2006 nearly filled the slit; as a re- 
sult the sky subtractions proved to be extremely difficult. 
Ultimately some of the star signals were removed along 
with the sky even when using the outermost 3-pixels on 
either side of the projected slit spectrum; thus, the final, 
co-added stellar spectra from 2006 (stars Pal l-I and -II) 
have lower SNR than expected. The sky subtraction for 
one of the 2007 stars, Pal 1-IV, proved impossible, and 
that star has therefore been dropped from further abun- 
dance analyses. The other two stars. Pal l-III and -C, 
have good sky subtraction and therefore higher SNR. 

The radial velocities were determined by locating sev- 
eral strong, easily-identified spectral lines and calculating 
the shift in wavelength using Gaussian fits to the lines. 
These radial velocities were then shifted to heliocentric 
velocities using the IRAF task rvcorrect; the final he- 
liocentric velocities are listed in Table [5] Though a sky 
subtraction could not be performed on Pal 1-IV, a ra- 
dial velocity could still be obtained; this radial velocity 
is quite different from the other stars, as will be discussed 
in Section FTOl 

The shifted rest-frame spectra were then continuum- 
normalized. The continuum was estimated using an it- 
erative, non-parametric filter (e.g. IBattaglia et al.|[2"008D 
designed to split the spectrum into continuum and line 
components based on a user-defined scale length for al- 
lowed continuum variations. This was set at 10 A so that 
the strongest spectral absorption lines were not affected. 
Figure |3] shows red and blue portions of the final stellar 
spectra which have a full coverage of approximately 4350 
to 7140 A; the representative spectra in Figure |3] have 
been arbitrarily shifted vertically for easier comparison. 

3. EQUIVALENT WIDTHS 

The measured s pectra l lines are from the line list 
by 'Shetrone et al.' (|2003|). with su p pleme nt al lines from 
.Cavrel et al.. (,2004i ). lAoki et al.l (|2007l ). iCohen et al.1 
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TABLE 1 

Our target stars, their positions, and their magnitudes 



Star 



RA (J2000) Dec (J2000) V^' 



r 



Kb 



References'^ 



M67-141 
Pal l-I 
Pal l-II 
Pal l-III 
Pal 1-IV 
Pal 1-C 



3'^ 33" 
3'^ 33" 
3'^ 33" 
31^33" 
3'^ 33" 



22.8= 
21.8= 
29.6= 
12.3= 
27.0= 
21.0= 



+11°48'1."7 
+79°35'16."2 
+79°34'16."2 
+79°34'59."2 
+79°35'34."9 
+79°34'57."1 



10.480 
16.675 
16.843 
17.827 
18.032 
16.603 



9.400 
15.459 
15.618 
16.628 
16.969 
15.328 



7.942 
13.832 
13.983 
15.281 
15.779 
13.715 



1,2, 
4 
4 
4 
4 
4 



^ M67-141 V and I magnitudes are in the Kron-Cousins system, wiiile the Pal 1 V and 
I magnitudes have been transformed to the Johnson system from the HST system. 

All K magnitudes are from the 2MASS Point Source Catalog. 
" References. (l)[H0g et al.l |[ 2000l) : (2) ISanderg 119771) : f3) IJanes fc SmithI fl98i) : (4) 
ISaraiedini et al.l 120071) 



TABLE 2 

Dates of observations, total exposure times, signal-to-noise ratios at two wavelengths, and 
calculated heliocentric radial velocities 



Star 



Dates 



Exposure Time (s) SNR^ (5310 A) SNR'' (6616 A) sjhcllo (km s"!) 



M67-141 


2006 Jan 21 




5985 




100 150 4-33.9 ±1.0 


Pal l-I 


2006 Jan 20, 21 




9000 




15 30 -77.2 ±1.0 


Pal l-II 


2006 Jan 20 




9000 




15 30 -78.0 ±1.0 


Pal l-III 


2007 Jan 25, 26 




12600 




31 45 -77.2 ±1.0 


Pal 1-IV 


2007 Jan 27, 28, 29 




21600 




-^ -b -68.1 ±1.0 


Pal 1-C 


2007 Jan 6, 27 




5400 




35 55 -77.0 ±1.0 


==■ SNR arc 


per resolution element. 










*= Pal 1-IV 


was dropped from the 


analysis. 


as the seeing was 


too 


poor to perform a sky subtraction. 



(l2008l). iLetarte et alj (|20Q9|) . iTafelmever eFal] (|20Tot ). 
and i Frebel et al.l (|2010D . The equivalent widths (EWs) 
were fit with Gaussian profiles using splot in IRAF and 
were checked by numerical interpretation under the con- 
tinuum. EWs can now be measured with automatic 
line-measuring programs such as DAOSPEClj a program 
that is intended for hi gh resolution (R > 1500 0), high 
SNR (> 30) spectra ( Stetson fc Pancinol [20081) . Some 
advantages of DAOSPEC over IRAF's splot include a 
fixed full-width at half maximum (FWHM) for all lines, 
and an effective continuum that takes weak features into 
account; however, it appears that with low SNR spec- 
tra DAOSPEC can have difficulty fitting the continuum. 
Due to the poor quality of the Pal l-I and -II spectra, 
IRAF splot measurements are therefore preferred. By 
a visual comparison with the higher SNR spectra (Pal 
l-III and -C), the spectral lines in Pal l-I and -II were 
identified and distinguished from noise in the low SNR 
spectra. Spectral lines in the noisy spectra that had a 
drastic difference in width or depth from Pal l-III or -C 
were discarded. 

Measurements were made relative to the local contin- 
uum. For all stars, lines stronger than 200 mA were 
rejected. For Fe, the Pal 1 lines with EWs stronger 
than 150 mA were also thrown out, in order to best 
constrain the atmospheric parameters (discussed in Sec- 
tion |3]). The wavelengths (in A), excitation potentials 
(in eV), log 5/ values, and equivalent widths (in mA) for 

5 DAOSPEC has been written by P.B. Stetson for the Domin- 
ion Astrophysical Observatory of the Herzberg Institute of Astro- 
physics, National Research Council, Canada. 



the four Pal 1 stars and the standard star are shown in 
Tabic H 

3.1. IRAF splot vs. DAOSPEC measurements 

DAOSPEC has been used for spectral line measure- 
ment s by several different authors (e.g. ILetarte et al.l 
l201Clf ): therefore we compare our splot EWs to those 
from DAOSPEC. Figure |4] shows a comparison between 
splot EWs (fro m this work) and DAOSPEC EWs (from 
IPancino et all [2010) for the standard star, M67-141, 
which has the highest resolution and SNR. Their EW er- 
rors from DAOSPEC are also shown. The majority of the 
offsets between the splot and the DAOSPEC measure- 
ments are within our adopted EW measurement errors 
(shown as dashed lines; see Section [5731 for a description 
of how these errors are calculated), and are therefore 
not significant. The few lines outside these errors are 
weak lines that we measure to be slightly stronger than 
DAOSPEC measures — this difference may be attribu ted 
to continuum placement fsee lStetson fc Pancinoll20()8l for 
a discussion of their effective continuum) . For the weak 
lines (> 60 mA), the average discrepancy between the 
two analyses (shown as a dotted blue line) is 4 ± 6 mA. 

3.2. A Comparison with Previous Studies 

Our EW measurements from the low SNR Pal l-I 
spectr um can be compared to those from [Monaco et al.l 
(|201Clf ). Unfortunately there are very few lines in com- 
mon for us to compare, due to differences in the spectral 
range and SNR . For the lines in common, Figure [5] shows 
our splot EWs versus the Monaco et al. DAOSPEC EWs 
from the higher SNR spectrum. Also shown are the line 
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TABLE 3 

The line list.^ 



Wavelength Element E.P. 
(A) (eV) 



log gf Equivalent width (niA) 

Sun M67-141 Pal l-I Pal l-II Pal l-III Pal 1-C 



4443.19 


Fel 


2.86 


-1.043 


134.0 


4476.02 


Fel 


2.85 


-0.819 


- 


4484.22 


Pel 


3.60 


-0.864 


100.0 


4489.75 


Pel 


0.12 


-3.899 


91.0 


4592.66 


Pel 


1.56 


-2.462 


99.0 



171.0 



122.0 



185.0 



76.0 
120.0 



157.0 
117.0 



^ Tabic [3] is published in its entirety in tiie electronic edition of The Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 




5570 5572 5574 

Wavelength (A) 
(a) Blue Region 




6018 6020 6022 6024 6026 6028 

Wavelength (A) 
(b) Red Region 

Pig. 3. — Portions of the five spectra, in both a blue and a red 
region, with notable spectral lines labelled. The spectra have been 
continuum normalized and were arbitrarily shifted vertically for the 
ease of comparison. Note the low S/N of Pal l-I and -II compared 
to the other Pal 1 stars. 



of perfect agreement (solid red line), our estimated la 
EW errors (black dashed lines; see Section [O)) . and the 
average trend (dotted blue line). Our EWs tend to be 
slightly higher, with an average offset of 7 ± 6 mA. How- 
ever, the majority of the points in Figure [5] lie within 
the dashed lines, which implies that the EW measure- 
ment techniques and spectral reduction methods are in 
reasonable agreement. We further note that there does 
not seem to be a significant trend in error with EW. 




20 40 



100 120 

EW (Pancino et al.) 



Fig. 4. — A comparison of the EWs as measured with 
in IRAP to those measured by the program DAOSPEC (from 
IPancino eFal. 2010 '! for the standard star M67-141. The red line 
shows equal EW measurements. The blue dotted line shows the 
average offset between the weak lines in the two data sets. The 
dashed lines s how the EW errors (EWmin + 10% EW) as calcu- 
lated using the lCavrel {1988) formula (see Section l5.3ll . 



4. MODEL ATMOSPHERES 

Stellar atmospheric parameters (effective temperature. 
Tee', surface gravity in its logarithmic form, log 5; micro- 
turbulence, ^, in km s~^; and metallicity, [Fe I/H]3) are 
initially derived from broadband photometry and revised 
based on spectral indicators as described below. 



4.1. Photometric Parameters 

Determining the atmospheric parameters from broad- 
band colors requires knowledge of the metallicity, redden- 
ing, distance modulus, and mass of the stars. The ini- 
tial metalli c ity for M67-141, [Fe/H] = 0.01, comes from 
lYong et"aLl ()2005[ ). while the initial metallicity for the 
Pal 1 stars, [Fe/H] = —0.6, come s from the CaT mea- 
surements of lRosenberg et aD (jl998b). The surface tem- 
peratures are estimated using the colors of the stars (from 
the magnitudes shown in Table[II. The ( V — K ) color to 
temperature conve rsion of lAlonso et al.l (|1999l with the 
lAlonso et al.ll200"ll correction) is used — this conversion is 

^ We use the standard notation: [X/H] = (12 -|- log(Nx/NH)) - 
(12 + log(Nx,0/NH,Q)) = loge(X) - log£(X)Q, where X is any 
element, and Nx and Nu are the column densities of element X 
and of hydrogen, respectively. 



Sakari et al. 



j.ovy 










• 1 ^_, ^- 
























•* ■■' / 












* ■"' /"^ 


160 










/^y' 














^ 120 
























^ 












■^ 100 










y.^ y 


^ 




















,* 


• y ^ 


60 






• _, 


> 


/• ^* 




(i) 


/. 


**•"" / 


^ 


^ 


40 






^ 
^ 







80 100 120 

EW (Monaco et al.) 



160 180 



Fig. 5. — A comparison betwee n our EWs for Pal l-I, and those 
measured bv lMonaco et afl I I2010I ) , for the li nes in common betwee n 
our analyses. EW errors were not given in [Monaco et all 1120101 ). 
The circled point is a Cr I line that is unusally discrepant; for 
that reason we have removed it from the analysis. T he d ashed 
lines represent the EW errors, as discussed in Section 15.31 The 
blue dotted line shows the average offset between the two data 
sets. Note that though we tend to measure slightly larger EWs, 
the differences are typically less than the EW errors. 



based on a large sample of Galactic GCs and field stars 
and is calibrated in the TCS photometric system. Con- 
versions fr om (F — jr)joh nson to {V — K')-YC?, are therefore 



used (Alonso et al. 



The lAlonso et al 



1991 ■ 
IT999I) 



calibrations are based on the 
absolute colors of the stars, meaning the cluster redden- 
ing must be taken into account. The adopted reddenings 
for Pal 1 and M67-141 are E{B - V) = 0.20 ± 0.03 and 
E{B-V) = 0.033±0.005. from the fSchleeel et al] ([l998h 
dust extinction maps. The E( V — K) r eddening is deter- 
mined with the conversions of lMcCaTll (|20b4). The final 
photometric temperatures are shown in Table |4l 

The bolometric correction for each st a r is calcu- 
lated using the formula from lAlonso et al.l ()1999[) , and 
the absolute bolometric magnitude is determined using 
the distance modulus. For M67-141 we si mply adopt 
the a bsolute distance modulus quoted in lYong et alJ 
(|2005f ): (m - M)o = 9.56. For Pal 1 there are 
two different literature valu es for the distance modu- 
lus: iRosenberg erall (Il9 98a1) give (m-M)o = 15.25 ±0.25 
whilelSarajcdinieiaLJlloO?) give (m-M)o = 15.65±0.10. 
Both are determined through main sequence fitting of 
Pal 1 to the lower main sequence of 47 Tuc, a Galac- 
tic globular c luster with a siin i lar m etallicity to Pal iH 
We adopt the lSarajedini et all ()2007l ) distance modulus, 
(m-M)y = 16.260 which is determined from HST pho- 
tometry. 

Finally, the photometric surface gravity is determined 
through a luminosity comparison with the Sun, which 
requires knowledge of the mass of our stars. All the 
targets are evolved stars: M67-141 is in the red clump 



'' As discussed in IRosenberg et al"] l ll998al ) , the presence of the 
sparse horizontal branch is uncertain, preventing accurate determi- 
nations of the distance modulus using zero age horizontal branch 
methods (e.g. IVandenBerg et al.ll2000fl . 

* We assume here that Ay = 3.07E{B - V) l|McCal]||200il . 



(|Yong et al.ll2005D and the Pal 1 stars are all presumably 
on the RGB. If we assume that there has been mini- 
mal mass loss since the stars evolved off the main se- 
quence, and if we assume that not much time has passed 
since the stars left the main sequence, then the main 
sequence turn off mass can be used as the RGB stellar 
mass. For most Galactic GCs this turnoff mass is typ- 
ically 0.8 M0 — however, as discussed earlier, both M67 
(an open cluster) and Pal 1 are younger than the typi- 
cal GC, and should therefore have a higher turnoff mass 
than 0.8 Mq. To determine the turnoff mass, we exam- 
ined t he Pal 1 isochrones computed by iSaraiedini et al.1 
(|2007f ): an age of 5± 1 Gyr gives a turnoff mass of Mto = 
1.14 zfc 0.06 M p). This estimate agrees with that from 
iCohenI ()2004D for Pal 12, another young globular cluster. 
M67-141 has an age from 4- 5 Gyr from isochrone fits 
(jVandenBerg fc Stetsonl '20041. which means the turnoff 
mass will be similar. The final photometric gravities are 
also shown in Table HI 



4.2. Spectroscopic Parameters 

Spectroscopic indicators are examined to refine these 
photometric parameters. The abundan ces of the ind i- 
vidual Fe lines are analyzed in MOOG (|Sneden|[l973l )FI 
using LTE OS MARCS model atmospheres with spheri- 
cal g eometries (jGustafsson et al.l 120081 IPlez fc Lambert] 
I2002f ). The final temperature was found by forcing the 
Fe I line abundances to be independent of the excitation 
potential (x); similarly, the microturbulence was refined 
by forcing the Fe I abundances to be independent of the 
line equivalent width. The associated errors in these pa- 
rameters were found using the la errors of the slopes of 
the Fe I abundances versus x and EW. 

It is not clear whether spectroscopic methods are ap- 
propriate for determining the gravity (i.e. ensuring that 
the Fe I and Fe II abundances are equal), since Fe I is 
expected to suffer from non-LTE effects which depopu- 
late this under represented species through interactions 
with the radiation field. In addition, there are few Fe II 
lines in our spectral range; those that are available are 
located at blue wavelengths so that reliable Fe II abun- 
dance determinations can be difficult, especially in low 
SNR spectra. For these reasons, we choose to adopt the 
photometric surface gravities and their associated errors. 

The final metallicity is the average Fe I abundance, 
with an error determined from the line-to-line scatter, 
cr(Fel), divided by the square root of the number of 

fines ((5(Fe I) = cr(Fe 1)/VN; see Section[5X2l) . The final 
spectroscopic values (Tgff, logg, £_, and [Fe/H]) for M67- 
141 and the Pal 1 stars are listed in Table ID The spec- 
troscopic and photometric temperatures are in excellent 
agreement. We also list the atmosp heric parameters from 
the literature for M 67-141 (from lYong et al.l 120051 and 
iPancino et al.1120101 ) and for Pal l-I (from [Monaco et al.l 

^ The 2010 version of MOOG was obtained at 
http://www.as.utexas.edu/~chris/inoog.html We note that 
the scattering version (MOOG-yCAT) was not used in this 
analysis because of the high metallicities of our targets. A test 
of MOOG-SCAT showed that our line abundances would have 
negligible corrections; the average correction for Pal l-III is <0.01 
dex for lines below 5000 A. 
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TABLE 4 
Atmospheric Parameters for the Pal 1 stars and the standard star, M67-141. 





Photometric 


Spectroscopic 








Star 


Teff (K) 


Teff (K) 


logg 


5 (km s-i) 


[Fe/H] 


vsini (km s ^)'^ 


M67-141: This paper 


4624 ± 12 


4750 ± 50 


2.31 ±0.06 


1.4 ±0.1 


0.02 ± 0.02 


4.5 


Yone et al. (20051 


4700 (4604) b 


4700 


2.3 


1.34 


0.00 


- 


Pancino et al. (2010) 


4590 


4650 


2.8 


1.3 


0.06 


- 


Monaco et al: Pal l-I 


4850 


5000 


2.40 


1.0 


-0.5 


_ 


Pal l-I 


4742 ±81 


4800 ± 70 


2.27 ±0.15 


1.4 ±0.3 


-0.61 ±0.08 


2.5 


Pal l-II 


4725 ± 80 


4750 ± 135 


2.33 ±0.15 


1.6 ±0.2 


-0.61 ±0.08 


2.5 


Pal l-III 


5061 ± 90 


5050 ± 50 


2.89 ±0.15 


0.9 ±0.1 


-0.60 ±0.02 


0.5 


Pal 1-C 


4698 ± 79 


4750 ± 70 


2.22 ±0.15 


1.3 ±0.1 


-0.58 ±0.03 


1.5 



^ This parameter was determined for spectrum sytheses, as discussed in Section |5.2| 

The temperature in parentheses is the photometric temperature with the spectroscopic metalhcity. 
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Fig. 6. — The abundance trends with excitation potential (top), 
reduced EW (middle) and wavelength (bottom) for Pal l-I, with 
the spectroscopic atmospheric parameters listed in Table |4] The 
black dashed line shows a constant trend, while the red dot-dashed 
line shows the trends of the data; note that in all three panels 
the trend is essentially constant, though there is still a significant 
dispersion. 



l201Clf ). O ur values for M67- 141 are in excellent agree- 
ment with lYong et al j (I2005D . but o ur temperature and 
gravity disagree with iPancino et aP (^2010). Our Pal l-I 
temperatures are slightly low er and our microturb ulence 
values are slightly higher than [Monaco et all (|2010[ ). The 
differences are primarily due to differing atomic data; as 
we explain in Section [5. 6. 21 these atmospheric parameter 
differences result in slightly different abundances in the 
various analyses. 

Figure [5] shows the abundances from the thirty Fe I 
lines that were measured in Pal l-I, generated using the 
final spectroscopic atmospheric parameters. Note that 
though there is little to no trend in abundance for excita- 
tion potential, reduced equivalent width, or wavelength, 
there is still a range of ~ 0.2 dex. 



5. ELEMENTAL ABUNDANCES 

The elemental abundances were determined in MOOG 
using EW matching. With these abundances, the loga- 
rithmic ratios with respect to Solar values, [X/Fe], are 



calculated. Tables [5] and H] show the [X/Fe] ratios, the 
statistical errors (see Section [?751 for a description of how 
these errors are determined), the number of lines used to 
determine the abundances, and comparisons with liter- 
ature values. These results will be discussed in Section 



5.1. Hyperfine Structure 



and 



The elements Sc, V, Mn, Co, Cu, Ba, La 
Eu were checked for hyperfine structure (HFS) correc- 
tions. Th e majority of the HFS and isotopic data 
come s from iProchaska et all (|2000l for Sc, V, Mn, and 



Co) , iBooth et al.l (119831. additio nal MnV iBiehl (Il976 , 
Cu),pcWilliam (1998, Ba), and iLawler et al.l (|2001al lb . 
La and Eu), with extra lines added from the Kurucz 
database^ The corrections per line were averaged over 
the entire star and were applied to the average abun- 
dances from MOOG. If HFS data were not available, 
then the line was neglected in the correction calculation. 
When the corrections were less than 0.1 dex they were 
ignored. Negligible corrections were found (< 0.1 dex) 
for Sc, Ba, La, and Eu; moderate corrections (< 0.5 dex) 
for V and Mn; and large corrections (< 1.0 dex) for Co 
and Cu. 

5.2. Spectrum Synthesis 

When there was only a single absorption line for an 
important element, a spectrum synthesis was performed 
on the line. In each spectrum synthesis the instrumen- 
tal and rotational broadening are taken into account. 
With resolutions of i? = 36000 and R = 60000, Pal 1 
and M67-141 have instrumental broadening values of 8.3 
and 5.0 km s~^, respectively. The rotational broadenings 
for each star were determined by examining nearby lines 
with known abundances, and adjusting the usini values 
until the widths and depths matched the points. For ex- 
ample, in the case of the the Eu II 6643 A line, the 6644 
A Ni I line and the 6647 A Fe I line were used to de- 
termine these parameters. The derived wsini values for 
M67-141 and the four Pal 1 stars are shown in Tabled 
The abundances were then estimated from the best fits 
to the adopted models, with an error in log e equal to the 
range of abundances that fit the line profiles. Where the 
spectral lines were particularly noisy, only upper limits 

^^ |http : //kurucz . harvard . edu/llnelists . html | 
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TABLE 5 

A COMPARISON OF THE DERIVED ABUNDANCES IN THE STANDARD STAR, 



M67-141 FROM THIS WORK. 


IYONGETAL.N2005I1. 


ANDIPANCINO et AL 


.1 


(120101) 


. We ALSO 


GIVE THE TOTAL ERROR IN [X/Fe] AND THE NUMBER OF 






LINES USED IN THE CALCULATION. 








Pancino et al. (2010') 


Yone et al. ('2005) 


This work 

[X/Fe]- 


S 




X 


[X/Fe]- 


s 


N 


Fel 


0.06 


0.01 


-0.01 


0.12 


0.02 


0.02 


91 


Fell 


0.01 


0.03 


0.01 


0.09 


-0.03 


0.04 


24 


OI 


-0.05 


0.09 


0.10 


0.06 


-0.15 (0.03) 


0.09 


2 


Nal 


0.10 


0.02 


0.24 


0.10 


0.29 


0.13 


4 


Mgl 


0.29 


0.03 


0.18 


0.04 


0.04 


0.07 


3 


All 


0.06 


0.06 


0.16 


0.06 


0.14 


0.15 


2 


Si I 


0.09 


0.02 


0.11 


0.08 


0.13 


0.05 


10 


Cal 


-0.13 


0.02 


0.09 


0.03 


0.07 


0.06 


15 


Sell 


-0.02 


0.08 


- 


- 


0.09 


0.06 


12 


Til 


-0.07 


0.02 


0.05 


0.05 


0.02 


0.05 


35 


Till 


-0.07 


0.02 


- 


- 


0.10 


0.05 


33 


VI 


0.13 


0.04 


- 


- 


-0.05 


0.05 


17 


CrI 


0.01 


0.03 


- 


- 


-0.04 


0.09 


8 


cm 


- 


- 


- 


- 


0.18 


0.18 


2 


Mnl 


- 


- 


-0.20 


0.03 


-0.08 


0.13 


4 


Col 


0.11 


0.02 


0.01 


0.09 


-0.15 


0.21 


2 


Nil 


0.06 


0.05 


0.06 


0.09 


0.08 


0.04 


26 


Cul 


- 


- 


- 


- 


0.16 


0.13 


3 


Znl 


- 


- 


- 


- 


0.05 


0.12 


3 


YII 


-0.04 


0.02 


- 


- 


-0.04 


0.11 


5 


Ball 


0.26 


0.05 


0.02 


- 


0.0 


0.17 


2 


La II 


0.06 


0.05 


0.13 


0.04 


0.04 


0.06 


4 


Ndll 


0.01 


0.29 


- 


- 


0.04 


0.08 


5 


EuII 


- 


- 


0.05 


- 


0.0 


0.12 


1 


^ [X/H] is given instead for Fe I 


and Fc ll. 











TABLE 6 

[X/Fe] ratios, 15 ERRO RS, AND THE NUMBER OF LINES USED, FOR THE PAL 1 STARS. THE [X/Fe] RATIOS FROM 
IMONACO ET AL.I ()2010l) ARE ALSO SHOWN FOR COMPARISON. 



Monaco Pal l-I 
[X/Fe]'' 



Pal l-I 

[X/Fe]^ 



N 



Pal l-II 

[X/Fe]'' 



N 



Pal l-III 

[X/Fe]'' <5 



N 



Pal 1-C 

[X/Fc]" 



N 



Fel 


-0.49 


-0.61 


0.08 


30 


-0.61 


0.08 


21 


-0.60 


0.03 


81 


-0.58 


0.03 


60 


Fell 


-0.53 


-0.46 


0.30 


1 


-0.64 


0.44 


2 


-0.61 


0.05 


23 


-0.68 


0.05 


23 


Ol 


- 


<0.82 


- 


1 


<0.32 


- 


1 


<0.42 


- 


1 


0.20 (0.23) 


0.11 


2 


NaT 


0.38 


0.26 


0.34 


2 


0.23 


0.31 


2 


0.20 


0.09 


4 


0.16 


0.09 


4 


Mgl 


0.11 


-0.11 


0.20 


1 


-0.13 


0.30 


1 


-0.02 


0.08 


5 


0.02 


0.07 


6 


All 


0.25 


- 


- 


- 


- 


- 


- 


0.0 


0.10 


2 


0.07 


0.10 


2 


Six 


-0.01 


0.24 


0.24 


1 


0.13 


0.23 


1 


0.09 


0.06 


8 


0.19 


0.06 


8 


Gal 


0.04 


0.16 


0.16 


6 


-0.04 


0.22 


4 


0.15 


0.06 


17 


0.12 


0.06 


18 


Sell 


-0.01 


0.28 


0.26 


3 


0.21 


0.33 


2 


0.30 


0.09 


8 


0.21 


0.07 


12 


Til 


0.10 


-0.03 


0.34 


3 


-0.14 


0.38 


4 


0.02 


0.06 


24 


-0.08 


0.06 


24 


Till 


- 


- 


- 


- 


-0.13 


0.51 


2 


-0.08 


0.05 


32 


0.08 


0.07 


23 


VI 


-0.06 


-0.05 


0.21 


1 


0.08 


0.14 


4 


0.14 


0.07 


6 


0.06 


0.05 


14 


On 


-0.23 


- 


- 


- 


- 


- 


- 


0.01 


0.09 


11 


-0.16 


0.11 


9 


Cm 


- 


- 


- 


- 


- 


- 


- 


0.17 


0.20 


2 


0.19 


0.21 


2 


Mnl 


-0.22 


- 


- 


- 


-0.16 


0.36 


1 


-0.21 


0.12 


4 


-0.11 


0.10 


8 


Col 


-0.06 


- 


- 


- 


- 


- 


- 


-0.10 


0.18 


2 


0.03 


0.20 


2 


Nil 


-0.03 


0.05 


0.25 


3 


0.09 


0.14 


5 


0.08 


0.05 


20 


0.03 


0.05 


21 


Cul 


- 


- 


- 


- 


-0.05 


0.33 


1 


0.04 


0.18 


2 


-0.08 


0.32 


1 


Znl 


0.38 


- 


- 


- 


- 


- 


- 


0.08 


0.14 


3 


0.15 


0.15 


3 


YII 


-0.32 


< 0.04 


- 


1 


<0.34 


- 


1 


-0.45 


0.12 


4 


-0.36 


0.12 


5 


Ball 


0.24 


0.27 


0.24 


2 


0.19 


0.26 


2 


0.22 


0.12 


3 


0.26 


0.18 


2 


La II 


0.29 


<0.38 


- 


1 


<0.78 


- 


1 


0.42 


0.16 


1 


0.24 


0.14 


3 


Ndll 


- 


0.19 


0.65 


1 


- 


- 


- 


0.13 


0.16 


2 


0.12 


0.14 


2 


EuII 


- 


<0.78 


- 


1 


<0.78 


- 


1 


0.50 


0.20 


1 


0.50 


0.20 


1 



[X/H] is given instead for Fe i and Fe ii. 
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5527 5528 5529 5530 

Wavelength (A) 

Fig. 7. — Spectrum synthesis for the Mg I 5528 A hne in 
Pal l-I Syntheses for three different Mg abundances are shown: 
loge(Mg) = 7.27, 6.97, and 6 .67. The parameters for the synthesis 
are discussed in Section 15.21 Note that though the fits match the 
width of the line well, the bottom-most point is not fit; this point 
is most likely a noise spike. 

are provided. Spectrum syntheses were also used when 
a known blend affected an important line, as in the case 
of the 6141.73 A Ba II li ne, which is blended with a Ni I 
line at solar metallicity (jAllende Prieto et al.ll200TI ). By 
using a spectrum synthesis the strength of the Ni line 
could also be taken into account. 

An example of a spectrum synthesis of the Mg I 5528 A 
line in Pal l-I is shown in Figure [T] Note that the width 
of the line is fit well, but the bottom point is missed; we 
suspect that this bottom point is a noise spike. 

5.3. Errors 

To determine the errors on the abundances, we sep- 
arate systematic errors (due to model atmosphere pa- 
rameters) from statistical errors (due to measurements). 
While the systematic errors are calculated and discussed, 
they are not included in the errors bars in our figures and 
are tabulated only in Tables [7] and HI 

5.3.1. Systematic Errors 

The abundances from individual spectral lines in a star 
are all affected by the choice of model atmosphere param- 
eters. Changing each parameter (Tofi, logg, ^, [M/H]) 
will affect the average abundance of each species. While 
the errors in the atmospheric parameters themselves are 
random, changes in these parameters will systematically 
affect the abundance analyses. 

To determine how the abundances were affected by un- 
certainties in the atmospheric parameters, new model at- 
mospheres were tested by varying one parameter at a 
time by its Icr error. These are quoted individually in 
Tables [7] and HI It is not clear how to combine these er- 
rors, since they are not independent, e.g. increasing \ogg 
affects the Toff and ^ estimates in our analysis. However, 
investigating their effects separately provides a maximum 
uncertainty estimate. We also consider systematic errors 
due to continuum placement; as discussed in Section [3l 
continuum placement can result in EW differences that 
are outside our adopted Icr errors. As a conservative 
estimate of these errors, we adopt the ave rage offset be 



(|2010l 7 mA, for the low SNR stars), for our continuum 
uncertainty. 

Representative systematic errors in abundances for 
M67-141, Pal l-I, and Pal l-III are shown in Tables [7] 
and m individually and added in quadrature. The dif- 
ferences are largest for the continuum and temperature 
errors in M67-141, while the errors are fairly similar for 
all parameters in the Pal 1 stars. The final errors in the 
[X/Fe] ratios are fairly small for all elements. 



5.3.2. Statistical Errors 



tween our EW measurements and those of Pancin o et al 
(|2blCl 4 mA, for the high SNR stars) and .Monaco et al 



Foll owing the method outlined by iShetrone et al.l 
(J20Q3!) , we identify three types of statistical errors for 
the abundances: the error from the Fe I line-to-line scat- 
ter, cr(Fe); the error from the line-to-line scatter of the 
element itself, (j{)^)] and the error in the abundances 
due to uncertainty in the measured equivalent widths, 
it(EW). Because there are quite a few Fe I lines, many 
of them strong enough to be detectable in the low SNR 
spectra, cr(Fe) should provide a conservative minimum 
estimate of the line-to-line scatter due to the SNR and 
continuum placement. If cr(X) is less than (j(Fe), we as- 
sume that (T has been underestimated because the spec- 
trum has too few lines measurable lines for element X. 
To esti mate t he errors due to EW measurements, the 
Cavrell (|1988f ) formula is used to determine the mini- 
mum EW measurement error (AEWmin) given the SNR 
and resolution per spectrum; the EW error for a single 
line is then A(EW)^ A(EW,„i„) + 10% EW (also see 
IShetrone et al.l 120031) . These EW errors are propagated 
through the abundance analyses, providing an average 
o'(EW). For each element, the larger of these three errors 
(cr(Fe), (t{X), cr(EW)) is selected, and divided by \/N\i-ncs 
(where iVnncs is the number of lines of that element) to 
determine the total statistical error, 5 (i.e. 5 = a/vN). 
The total 16 errors in [Fe/H] and [X/Fe] are listed in 
Tables [5] and El 

5.4. Solar Abundances 

An analysis was also performed of the solar spectrum 
from the Kurucz 2005 solar flux atlasE3 The model at- 
mosphere was generated with Tcs = 5777 K, logq = AAA 
dex, £ = 0.85 km s-\ and [M/H] = 0.0 (jYong et al.l 
120051) . Lines were discarded if they were too weak, too 
strong, or looked to be an unresolved blend at the solar 
Toff and [Fe/H] . The fin al solar abundances a re shown in 
Table M along with the lAsplund et all ()2009f ) values for 
reference. The results are in excellent agreement, though 
O, Sc, and Ba are all higher in our anal ysis, a trend that 
was also noticed by lYong et al.l (120051 for O and Ba) , 
IPancino et all ()20ia for Ba) and IMonaco et al.l (|2010l 
for Sc and Ba). By using these high solar values, sys- 
tematically high [0/Fe], [Sc/Fe], and [Ba/Fe] values will 
be removed from our Pal 1 analysis. Because the single 
Eu line is ver y weak in the solar spectrum, we simply 
use the Asplu nd et al.l ([2009) value. Our solar values are 
adopted in all calculations of [Fe/H] and [X/Fe] . 



|http : //kurucz . harvard . edu/sun . html] 
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TABLE 7 
Errors in M67-141's [X/Fe] ratios due to atmospheric parameter errors. The total error shows the 

individual errors added in quadrature. 



AT = +50 K Alog3 = +0.06 A^ = +0.1 km s" 



A[M/H] = +0.02 Continuum +4 mA Total 



Fc I/H] 


+0.07 


+0.01 


-0.05 


Fc II/H] 


-0.07 


+0.03 


-0.05 


O I/Fe] 


-0.06 


+0.02 


+0.05 


Na I/Fe] 


+0.01 


-0.01 


+0.02 


Mg I/Fe] 


-0.03 


0.0 


+0.04 


Al I/Fe] 


0.0 


-0.01 


+0.03 


Si I/Fe] 


-0.10 


0.0 


+0.03 


Ca I/Fe] 


+0.04 


-0.01 


0.0 


Sc Il/Fe] 


-0.08 


+0.01 


0.0 


Ti I/Fe] 


+0.09 


-0.01 


-0.01 


Ti Il/Fe] 


-0.08 


+0.02 


-0.01 


V I/Fe] 


+0.11 


0.0 


+0.01 


Cr I/Fe] 


+0.08 


-0.01 


-0.02 


Cr Il/Fc 


-0.13 


+0.01 


-0.01 


Mn I/Fc 


+0.05 


0.0 


-0.02 


Co I/Fc] 


0.0 


0.0 


-0.02 


Ni I/Fe] 


-0.03 


0.0 


0.0 


Cu I/Fe 


0.0 


-0.01 


-0.01 


Zn I/Fe 


-0.11 


+0.01 


0.0 


Y Il/Fc 


-0.07 


+0.02 


-0.02 


Ba Il/Fc 


-0.04 


+0.02 


-0.03 


La Il/Fc 


-0.05 


+0.02 


+0.04 


Nd Il/Fe 


-0.06 


+0.02 


+0.02 


Eu Il/Fc 


-0.08 


+0.02 


+0.04 



-0.01 


+0.07 


0.11 


-0.01 


+0.09 


0.13 


0.0 


+0.02 


0.08 


0.01 


-0.01 


0.03 


0.0 


0.0 


0.05 


0.01 


-0.01 


0.03 


0.01 


0.0 


0.10 


0.01 


-0.01 


0.04 


0.0 


0.0 


0.08 


0.01 


0.0 


0.09 


0.0 


+0.02 


0.09 


0.0 


+0.01 


0.11 


0.01 


0.01 


0.08 


0.01 


+0.02 


0.13 


0.0 


0.0 


0.05 


0.01 


+0.01 


0.02 


0.01 


0.0 


0.03 


0.0 


-0.02 


0.02 


0.0 


+0.03 


0.11 


0.0 


+0.03 


0.08 


-0.01 


0.0 


0.05 


0.0 


+0.02 


0.07 


0.0 


+0.05 


0.08 


0.0 


+0.02 


0.09 



TABLE 8 
Errors in Pal l-I and -Ill's [X/Fe] ratios due to the atmospheric parameter errors. The total error 

SHOWS THE INDIVIDUAL ERRORS ADDED IN QUADRATURE. 





ATeH 


(K) 


A log 3 


A^ (km s-i) 


A [M/H] 


Continuum (mA) 


Total 




I 


HI 


I 


HI 


I 


HI 


I 


HI 


I 


HI 


I 


HI 




+70 


+50 


+0.15 


+0.15 


+0.3 


+0.1 


+0.08 


+0.02 


+7 


+4 






A 


Fe I/H] 


+0.07 


+0.05 


0.0 


-0.01 


-0.18 


-0.04 


+0.01 


0.0 


+0.12 


+0.08 


0.23 


0.10 


A 


Fc II/H] 


-0.04 


-0.03 


+0.05 


+0.06 


-0.08 


-0.04 


+0.02 


0.0 


+0.11 


+0.14 


0.15 


0.16 


A 


O I/Fe] 


- 


-0.04 


- 


+0.08 


- 


+0.04 


- 


+0.01 


- 


+0.14 


- 


0.17 


A 


Na I/Fe] 


-0.01 


-0.02 


-0.02 


-0.01 


+0.08 


+0.02 


-0.01 


0.0 


-0.01 


-0.03 


0.08 


0.04 


A 


Mg I/Fe] 


-0.02 


-0.02 


-0.05 


-0.02 


+0.08 


+0.03 


0.0 


0.0 


-0.06 


-0.02 


0.11 


0.05 


A 


Al I/Fe] 


- 


-0.02 


- 


+0.01 


- 


+0.04 


- 


0.01 


- 


+0.03 


- 


0.05 


A 


Si I/Fe] 


-0.08 


-0.04 


+0.02 


+0.03 


+0.12 


+0.04 


0.0 


+0.01 


-0.01 


+0.02 


0.15 


0.07 


A 


Ca I/Fe] 


0.0 


0.0 


-0.02 


-0.02 


+0.03 


+0.01 


-0.01 


0.0 


-0.02 


-0.04 


0.04 


0.05 


A 


Sc Il/Fc] 


-0.08 


-0.05 


+0.06 


+0.07 


+0.09 


0.0 


+0.01 


0.0 


+0.01 


+0.05 


0.14 


0.10 


A 


Ti I/Fc] 


+0.04 


+0.02 


+0.01 


+0.01 


0.0 


0.0 


-0.01 


0.0 


+0.04 


+0.02 


0.06 


0.03 


A 


Ti Il/Fc] 


- 


-0.05 


- 


+0.06 


- 


0.0 


- 


0.0 


- 


+0.05 


- 


0.09 


A 


V I/Fe] 


+0.05 


+0.02 


0.0 


+0.02 


+0.14 


+0.03 


-0.01 


0.0 


0.0 


+0.03 


0.15 


0.05 


A 


Cr I/Fc] 


- 


+0.02 


- 


-0.01 


- 


-0.01 


- 


0.0 


- 


-0.01 


- 


0.03 


A 


Cr Il/Fc 


- 


-0.07 


- 


+0.07 


- 


0.0 


- 


+0.01 


- 


+0.07 


- 


0.12 


A 


Mn I/Fc 


- 


+0.01 


- 


-0.01 


- 


0.0 


- 


+0.01 


- 


-0.02 


- 


0.03 


A 


Co I/Fc] 


- 


0.0 


- 


+0.03 


- 


+0.02 


- 


+0.01 


- 


+0.04 


- 


0.05 


A 


Ni I/Fc] 


-0.01 


-0.02 


+0.02 


+0.01 


+0.01 


0.0 


-0.01 


0.0 


+0.01 


+0.01 


0.03 


0.02 


A 


Cu I/Fc 


- 


0.0 


- 


+0.01 


- 


-0.02 


- 


0.0 


- 


+0.01 


- 


0.02 


A 


Zn I/Fc 


- 


-0.07 


- 


+0.04 


- 


0.0 


- 


0.0 


- 


+0.05 


- 


0.09 


A 


Y Il/Fe 


- 


-0.05 


- 


+0.07 


- 


0.0 


- 


+0.01 


- 


+0.08 


- 


0.12 


A 


Ba Il/Fe] 


-0.05 


-0.04 


+0.05 


+0.04 


-0.08 


-0.02 


+0.02 


+0.01 


0.0 


-0.01 


0.11 


0.06 


A 


La Il/Fe] 


- 


-0.04 


- 


+0.07 


- 


+0.03 


- 


0.0 


- 


+0.10 


- 


0.13 


A 


Nd Il/Fc 


-0.05 


-0.05 


+0.07 


+0.07 


+0.09 


+0.01 


+0.01 


0.0 


+0.05 


+0.08 


0.13 


0.12 


A 


Eu Il/Fe 


- 


-0.05 


- 


+0.08 


- 


+0.04 


- 


+0.01 


- 


+0.14 


- 


0.17 
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Fel 


7.45 


Fell 


7.45 


OI 


8.69 


Nal 


6.27 


Mgl 


7.53 


All 


6.43 


Si I 


7.51 


Cal 


6.29 


Sell 


3.05 


Til 


4.91 


Till 


4.91 


VI 


3.96 


CrI 


5.64 


CrII 


5.64 


Mnl 


5.48 


Col 


4.87 


Nil 


6.20 


Cul 


4.25 


Znl 


4.63 


YII 


2.17 


Ball 


2.18 


La II 


1.17 


Ndll 


1.45 


EuII 


0.51 



TABLE 9 

Our solar a bundances compared t o those 

from lasplund et al.i ii2009i) . 

Element Asplund et al. (2009) This work 

7.50 

7.49 

8.83 (8.85)'' 

6.25 

7.59 

6.47 

7.51 

6.32 

3.19 

4.93 

4.95 

3.95 

5.63 

5.59 

5.41*= 

4.91'' 

6.27 

4.28'' 

4.58 

2.18 

2.29 

1.22 

1.55 

0.5r 



^ The number in parentheses is the O abundance as 

determined with MOOG 2002. 

^ With HFS corrections 

'^ From Asplund et al. (2009) 

5.5. NLTE Effects 

Non-LTE corrections are neglected in this analysis, yet 
it is well known that many of our elements do suffer from 
NLTE effects. We have tried to minimize this problem 
by eliminating lines which are greatly affected by NLTE. 
Na is a particularly sensitive element; however, weaker 
lines tend to ha v e sma ller corrections. At solar metal- 
hcitv lLind eFaH (|201ir ) recommend the 6154.23/6160.75 
A doublet, with the 5682.65/5688.21 A doublet for 
[Fe/H] < -1.0. For the Sun and M67-141 only the for- 
mer doublet is used; for Pal 1, however, both sets are 
used. Any c orrections to these lines should be < —0.10 
dex at most ("Mash onkina et al1l2000( l. 

Th e NLTE effects for Mn c an als o be quite large. How- 
ever, iBergemann fc GehrenI (|2007| ) note that the correc- 
tions are strongest for the weak and intermediate lines 
(i.e. those with EW < 80 mA). All of our lines are 
stronger than 80 mA — such corrections are then "scat- 
tered around zero or negative." 

Barium is another element that suffers from NLTE ef- 
fects that can broaden the lines. The 6496.91 A line 
has been eliminated from this analysis, since it has large 
corrections. However the corrections for 5853.69 A and 
6141.73 A are negligible for Pal I's metallicity range, 
and the 4554.03 A correction is small at solar metallic- 
ity (IShort fc Hauschildtll2006h . 

5.6. Comparisons with Previous Work 

We compare our [X/Fe] results for M67-141 to 
lYong et all (I2005D andlPancino et all (|2010D . and for Pal 
l-I to [Monaco et al.l (|2010[ ). in Tables [5] and HI Here we 
discuss any discrepancies. 



5.6.1. M67-141: Yong et al. (2005), Pancmo et al. (2010) 

Our [X/Fe] va lues agree very well with those of 
lYong et all ()2005l see Table [SJ note that we do not shift 
their values to common solar abundances as they de- 
rive their own solar abundances). The only inconsis- 
tent abundances are those of O and Mg. The Mg dif- 
ference may be due to the choice of lines. The O abun- 
dance, however, relies on molecular equilibrium calcu- 
lations within MOOG; these calculations seem to differ 
between the 2002 version used by Yong et al. and the 
2010 version we use. Using the 2002 version increases our 
O abundances and brings the results into better agree- 
ment. These abundances are shown in Tables [5l[6l and 
[HI in parentheses following the 2010 versions. Note that 
the solar abundance is only slightly affected. As many of 
the comparison studies in Section |6] use the 2002 version 
of MOOG, we consider that version to be more accurate 
for comparisons. 

Our [X/Fe] ratios do not agree as well with 
iPancinoet alTllMol) : Na, Mg, Ca, Ti, V, Co, and Ba 
all disagree. For Na, Al, Ti, and Ba the discrepancies 
seem to be due primarily to the different atmospheric 
parameters. iPancino et al.l (|201Cl( ) also do not appear to 
have applied HFS corrections to V or Co, which could ex- 
plain why their abundances are much higher than ours. 
Mg is again discrepant because of the choice of lines and 
atomic data. Ca is also quite low in their analysis, also 
apparently due to atomic data. 

Finally, we note that a recent spectral comparison by 
lOnehag et al.l (|2010f ) of a solar twin in M67 suggests that 
the entire cluster should have roughly solar abundances, 
with a slightly elevated [Fe/H] ~ 0.02 dex. With the 
exceptions of Na and Si, all of our elements are within 
0.08 dex of solar [X/Fe]. Na c ould very well be e l evated 
in clump stars like M67-141; iTautvaisiene et al.l ()2000f ) 
suggested that mixing could take place in open clusters 
after the He-flash, which could bring up material from 
the Na-Ne cycle. Alternatively this high Na abundance 
could indicate a need for a negative NLTE correction. 



5.6.2. Pal l-I: Monaco et al. (2010) 

As shown in Table [HI Monaco et al.'s analysis differs 
from ours for Mg, Al, Sc, and Ti. The significantly dif- 
ferent atmospheric parameters affect most of these ra- 
tios, particularly Al (since our Pal l-I spectrum does 
not have a sufficiently high SNR to determine the Al 
abundance, we compare their Pal l-I abundance to our 
Pal l-III and -C values, which assumes there is not an 
abundance spread). In addition, different atomic data 
for Mg, Al, and Ti causes additional discrepancies. The 
Sc II abundance seems to differ due to the solar Sc abun- 
dances; the measured solar EWs in lMonaco et al.l ()2010f ) 
are considerably higher than ours, leading to a higher 
solar Sc abundance and therefore a lower [Sc/Fe] for Pal 
l-I. 

The differences in atomic data are rather discouraging; 
however, we note that while the choices of lines are broad, 
much of ou r atomic data overlaps wit h that of lFulbrighti 
(|2000l) and iReddv et al.) (|2003l . |2006| ). who provide the 
majority of the metal-rich Galactic stars for the com- 
parisons in our analysis. Therefore we consider that our 
results are fairly robust for this comparison. 
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6. ABUNDANCE RESULTS 

Figures [S] through [12] compare the abundance trends 
of the Pal 1 stars (red star symbols) and M67- 
141 (red circle) to those of Galactic field stars (field 
stars; grey poi nts, from the sour ces assembled by 



iVenii et al.N2004 with add itions from lReddvet all 120061 
and iSimmerer et alJ 120041 ) . Galactic GCs (black solid 
symb ols, fr om the sources assembled in Pritzl et alj 
[2001 plus ISbordoneetall [20071 and Cohen 20oC 
and old, metal-poo r Galactic P Cs (op en ci rcles, from 
Carraro et all 12004 lYong et all 120051 and iFriel et alJ 



2010| ). The two old, most metal-poor OCs are Berke- 



ley 20 (Be 20) and Be 29; the latter has been linked 
to the Sgr dSph based on it s location, chemistry, and 
kinematics (jCarraro &: Bensbv Il2009|) and its location o n 
the Sgr age-metallicity relation (jForbes fc Bridge J2010I ). 
The bulge and thin/thick disk GCs (circles) are distin- 
guished from the halo clusters Pal 12 (plus signs) and 
Ter 7 (crosses) , both of which are associated with the Sgr 
dSph. The bulge and disk cluster abundances are aver- 
aged over the entire cluster, while the individual stars in 
the halo clusters are separate. Those studies that do not 
perform their own solar abundanc e analyses have been 
shifted to the lAsplund et al] ()2009[ ) abundances. 

6.1. Metallicity 

The derived metallicities are in excellent agreement 
with previous studies. For M6 7-141, [Fe/H] = .02 ± 
0.02, which agrees with both lYong et al.l l|2005D and 
lOnehag et al.l (j2010|) . The average metallicity of the 
four Pal 1 stars, [Fe/H] = — 0.60 ± 0.01, is in agree - 
ment with its CaT metallicity ([Rosenberg et al.lll998bD 
and i ts metallicity fro m isochrone fit s (iSar aiedini et al.l 
|2007| ). As discussed bv iMonaco et al . (2010), this aver- 
age metallicity is high for a GC so far from the Galactic 
center — the only metal- rich Galactic halo GCs are Pal 12 
and Ter 7, neither of which are actually Galactic. How- 
ever, there are bulge and disk clusters in this metallicity 
regime, including 47 Tuc. This metallicity is also low for 
an OC, though Be 20 and 29 are nearly as metal-poor. 

We further note that the [Fe I/H] and [Fe II/H] ratios 
are in rough agreement, suggesting that our results would 
not differ if spectroscopic gravities had been used (see 
Section H]). 

6.2. a-Elements 

The a-elements form via the capture of '*He nuclei. 
This happens primarily during hydrostatic burning in 
massive stars, and the a-elements are subsequently re- 
leased into the interstellar medium by Type II super- 
novae. Figures |8| and |9| show the [X/Fe] ratios versus 
[Fe/H] for a-elements (O, Mg, Si, and Ca) and elements 
that behave like a-elements in Galactic halo stars (Ti). 
Ti is an average of Ti I and Ti II, except for Pal l-I, 
which includes only Ti I — this average is not weighted 
by the number of lines. The total a abundance shown 
in Figure |9| is an average of Mg, Ca, and Ti. O was ex- 
cluded from this average because of its weak lines; Si was 
excluded because our [Si/Fe] may be systematically too 
high (see Sectionl5.6p . Note that the [a/F e] values from 
iVenn et al.l (|2004[ ) and iPritzl eFall ([2005D have been re- 
calculated using our definition. 

The O I lines in our spectral range are the forbidden 
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Fig. 8. — [X/Fe] ratios of o-elements O, Mg, and Si versus metal- 
licity for the stars in our study as well as Galactic field stars and 
clusters. The red star symbols are the four Pal 1 stars. M67-141 is 
the red circle. Points with arrows represent upper (or lower) limits. 
The grey points ar e Gala ctic field stars, from the sources assem- 
bled bylVenn ot al] |[200l ) with additions from Rcddy et al. (2006). 
Black op en circles are Galactic OCs fro m Carraro et al. (2004), 
lYong et al., f^OS), and IFriel eFaLl (1201^) . Black filled circles are 
bulge and disk GCs from the sources assembled by j-'ritzl ct al] 
IJ2005I ). Plus signs are stars in Palo mar 12 ffrom [Cohenl 1200^ ) 
while crosses are stars in Ter 7 (from ISbordone et al.ll2007l) : both 
are halo GCs. The dashed lines show solar values. 
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Fig. 9. — [X/Fe] ratios of the o-elements Ca and Ti, as well as 
of the average a from Mg, Ca, and Ti. The Ti abundance is an 
average of Ti I and Ti II, except for Pal l-I, whose Ti I abundance 
is shown since its Ti II abundance is abnormally high. The points 
are as in Figure |8] 



6300 and 6363 A transitions. Both of these lines are ex- 
tremely weak and had to measured with spectrum syn- 
theses. Only upper limits are available for Pal l-I, -II, 
and -HI. Still, the single O abundance seems to follow the 
general trend of the other a-elements. Pal 1-C's slightly 
low [0/Fe] value agrees with the lower range of Galactic 
field star abundances and with the old OCs, particularly 
Be 20 and 29. 

The IVIg abundances were determined from lines of 
varying strengths for 1VI67-141, Pal l-III, and -C. 1VI67- 
141 agrees nicely with the field stars. Only one line was 
detectable in Pal l-I and -II, and had to be measured with 
spectrum syntheses. This 5528 A line is very strong in 
all the stars, though it is under the 200 mA limit for the 
Pal 1 stars. Given that the Pal l-I and -II abundances 
agree with -HI and -C, they seem reliable. Overall the 
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[Mg/Fe] ratios in Pal 1 are lower than the Galactic field 
stars, the bulge and disk GCs, and the old OCs, but in 
agreement with the halo clusters Pal 12 and Ter 7. 

The Si abundances were also determined from a wide 
range of lines, but again only a single line (at 6155 A) was 
detectable in Pal l-I and -II. This line is of intermediate 
strength in the other two Pal 1 stars, but again the Pal 1- 
I and -II abundances agree with the others. As discussed 
in Section [521 M67-141's [Si/Fe] ratio is slightly higher 
than previous values in the literature; this suggests that 
Pal I's [Si/Fe] might also be slightly high. Its current 
[Si/Fe] value places it in agreement with the Galactic 
field stars, the halo GCs, the bulge/disk cluster 47 Tuc, 
and the old OCs. 

The Ca abundance was determined from four to eigh- 
teen lines for the Pal 1 stars, with Pal l-II having the 
lowest number. All four of these Ca I lines are fairly 
strong (> 100 mA) in Pal l-II, which has a slightly lower 
Ca abundance than the other three stars. Again these 
[Ca/Fe] values agree with Galactic field stars, GCs, and 
OCs, and the Sgr GCs. M67-141, with 15 Ca hues (most 
of them strong) matches the [Ca/Fe] values found in the 
field stars. 

A wide assortment of Ti I and II lines were used, of 
varying strengths. Only three to four Ti I lines were de- 
tectable in the Pal l-I and -II spectra; those lines were 
of strong and intermediate strength. The other spectra 
had considerably more Ti I lines. Because many of the 
Ti II lines are located in the blue regions of our spectra, 
which have lower SNR, no Ti II lines are detectable in 
Pal l-I, and only two intermediate strength lines are de- 
tectable in Pal l-II. Even with only two lines. Pal l-IFs 
Ti I and II abundances are in excellent agreement, as 
are those of Pal l-III (with 32 Ti II lines), -C (with 23), 
and M67-141 (with 33). The Ti abundances for M67-141 
again follow the field star trend, though Pal I's [Ti/Fe] 
ratios are slightly lower than Galactic field stars, GCs, 
and OCs. The Pal 12 and Ter 7 [Ti/Fe] ratios are also 
slightly higher than Pal 1. 

Overall, M67-141's a abundance ratios agree well with 
the Galactic field stars. Pal 1 shows good agreement with 
the halo GCs Pal 12 and Ter 7. With the exception of 
Si and Ca, the [X/Fe] ratios for the a-elements in Pal 1 
are slightly lower than the bulge and disk clusters and 
the Galactic field stars. The Galactic [a/Fe] ratios are 
distinct between the thick disk and the thin disk stars; 
Pal 1 is clearly lower than the thick disk stars, but is 
not clearly separated from the thin disk. Thus, Pal I's 
[a/Fe] ratios are not distinct from the Galactic field stars, 
GCs, and OCs within the 15 errors^^ however, we note 
that the Pal 1 stars are generally lower than the average 
Galactic trend. 

6.3. Na and Al 

Na and Al are also produced in massive stars during 
nucleosynthes is through carbon burnin g and hydrogen 
sheU burning (jWooslev fc WeavedllQQSh . The [X/Fe] ra- 
tios are shown in Figure [TOl Note that Galactic GCs 
all show star-to-star variations in Na and Al, and the 
average abundances of the GCs may not reflect the ac- 
tual primordial abundances. As discussed in Section 

12 Recall that 15 = lcj/y/N^~ (see Section [53^ . 
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Fig. 10. — [X/Fe] ratios of the light elements Na and Al versus 
metallicity. The points arc as in Figure [8] 



15.51 the Na lines with large NLTE corrections have been 
avoided. Thus, M67-141's Na abundance is based on the 
6154/6161 A doublet; the Pal 1 analysis also includes 
the 5682/5688 A doublet. Both Pal 1 and M67-141 have 
slightly higher [Na/Fe] values than most of the field stars, 
in agreement with the old OCs and GCs. The Pal 1 
[Na/Fe] are distinctly higher than in the halo GCs Pal 
12 and Ter 7. 

The Al abundances come from the 6696/6698 A lines, 
both of which are weak in the Pal 1 spectra — thus, 
[Al/Fe] ratios are not found for Pal l-I and -II. M67-141's 
[Al/Fe] agrees with the field stars, while Pal I's [Al/Fe] 
is slightly low for field stars, bulge/disk GCs, and OCs. 
As with sodium above, the [Al/Fe] ratios in the Pal 1 
stars are distinct from those of stars in Ter 7. 

Anticorrelations between Na/ and Al/Mg a r e seen 
in nearly all Galactic GCs, and iCarretta et al.l (|2010l ) 
suggest that all bona fide GCs show a range in these 
elements. In Pal 1 Al is clearly not enhanced, and is 
slightly low for Galactic stars; without significant Al 
overabundances Mg depletion is not expected to occur 
(jCarretta et al.l[2009[ ). Na, however, is slightly enhanced 
in Pal 1, lying just above the Galactic field stars. In 
the absence of any stellar evolutionary effects, Na should 
behave like the a-elements, which are also distributed 
by Type II supernovae. Thus Pal I's high Na abun- 
dances may indicate that Na-Ne cycled gas is present — 
since Al is not enhanced, this suggests that Al-Mg cycled 
gas is not present. With only upper limits on [0/Fe] for 
three stars, we can only comment that Pal 1-C's [0/Fc] is 
slightly low, as expected with higher [Na/Fe] . We further 
note that there is not a significant range in the [Na/Fe] 
values for the four stars (see Section [7.1.ip which suggests 
that the high Na is not a sign of the canonical Na/0 an- 
ticorrelation seen in Galactic GCs — however, we cannot 
rule this out with only four stars. 

6.4. Iron-Peak Elements 

The iron-peak elements are formed in both Type la and 
Type II supernovae, though the precise contributions of 
the different types are unknown. In Pal I's metallicity 
regime T ype la supernov a e shou ld be the dominant con- 
tributors; llwamoto et al.l ()1999D estimate that a Type II 
supernova may create ~ 0.1 Mq of iron-peak material 
while a la may contribute ~ 0.8 Mq. Figures ITT] and [T2l 
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Fig. 11. — [X/Fe] ratios of the iron-peak elements Sc, V, and Mn 
versus metallicity. The points are as in Figure |8] 

show the [X/Fe] ratios of the hon-peak elements. Note 
that in the single plot for [Cr/Fe] both Cr I (red) and 
Cr II (yellow) are shown. 



6.4.1. Scandium, Vanadium, Manganese, and Cobalt 

The odd-Z elements are similar because they have only 
a single stable isotope each; these elements also require 
corrections for hyperfine splitting. Note that Sc is in- 
cluded in this discussion, even though it is not a tradi- 
tional iron-peak element, as its production site is primar- 
ily Type II supernovae. 

Sc II and V I have various detectable lines in our spec- 
tra. Spectrum syntheses were performed on the 5527 A 
Sc II line for Pal l-II and on the 6126 A V I line for -I. 
Three intermediate strength Sc II lines (5667, 5669, and 
6005 A) were detectable in Pal l-I, while four interme- 
diate strength V I lines (6090, 6150, 6199, and 6225 A) 
were found in Pal l-II. 

Only four Mn I lines were detectable and under the 200 
mA limit in M67-141, all of them strong (> 100 mA). As 
discussed in Section [531 any NLTE corrections to these 
lines should be small. No Mn I lines were measurable in 
Pal l-I, and only a single line (at 6022 A) was measured 
in Pal l-II. The other two Pal 1 stars all have Mn I lines 
with EWs < 80 mA, which should require small NLTE 
corrections. 

Co I has only two lines in this spectral range, at 5483 
and 5647 A. Neither line was detected in Pal l-I or -II. 
In M67-141, Pal l-III, and -C the 5483 A hue is in- 
termediate or strong while the 5647 A line is weak to 
intermediate. 

All of Pal I's odd iron-peak elements agree with the 
Galactic field stars. Sc is slightly higher than the general 
Galactic trend, putting Pal 1 in closer agreement with 47 
Tuc than with Ter 7 or Pal 12. The V abundances are 
slightly lower than the Galactic trend, in agreement with 
field stars and Ter 7. [Mn/Fe] is slightly higher than all 
the GCs, but agrees with field stars and the old OCs Be 
20 and 29. Finally, [Co/Fe] is slightly low, in agreement 
with the halo clusters but still within range for OCs and 
field stars — a comparison cannot be made with the GCs 
as Co abundances are not available for all clusters. 
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Fig. 12. — [X/Fe] ratios of the iron-peak elements Co, Cr, and 
Ni versus metallicity. Both Cr I and Cr II are shown (Cr I in red, 
Cr II in yellow); note that the low SNR stars Pal l-I and -II only 
have Cr I abundances. The points are as in Figure [8] 

6.4.2. Chromium and Nickel 

There are 11 available Cr I lines in our line list, but 
only two Cr II lines, both of which are in the blue spectral 
regions. Thus, Cr II abundances are not available for Pal 
l-I and -II; unfortunately, the redder Cr I lines were also 
not measurable in Pal l-I and -II. For the other stars the 
Cr I abundances are in good agreement with each other 
and with the literature. M67-141 looks like Galactic field 
stars and GCs, while Pal 1 lies slightly below the field 
stars and agrees with the halo GCs. However, the Cr II 
abundances are considerably higher (identified as yellow 
in Figure [T2)) , an effect that has been observed in previous 
studies of Galactic stars (e.g. Zhang et al. 2009). The 
Cr II abundances may not be reliable because the lines 
are in the blue regions of our spectra; however, the dis- 
crepancy bet\ Yeen Cr I and Cr II may be due to NLTE ef- 
fects for Cr I (jSobeck et al.ll2007t [Bergemann fc Cescuttil 
I2010D . Therefore, neither [Cr I/Fe] nor [Cr II/Fe] may be 
valuable choices for comparisons. 

Ni also has many spectral lines, though only three and 
five are detectable in Pal l-I and -II, respectively. Both 
M67-141 and Pal I's [Ni/Fe] values match those of field 
stars and clusters, with the exceptions of Ter 7 and Pal 
12, which are slightly lower. 

6.5. Copper and Zinc 

The nucleosynthetic origins of Cu and Zn are varied, 
and the precise yields from various sites depend on the 
models used. Explosive burning in Type la and II su- 
pernovae can create Cu and Zn, as can neutron cap- 
tures. Given the uncertainty surrounding the formation 
of these elements they are discussed separately from the 
rest. Their [X/Fe] ratios are shown in Figure [T^ 

Three Cu I lines are examined here: 5105.5, 5700, and 
5782 A. These lines are all strong in M67-141, though 
the last two are < 100 A in the Pal 1 stars. None are 
detectable in Pal l-I, while only the last is detectable in 
Pal l-II. The 5782 A line was just off the end of the red 
region in Pal l-III and Pal 1-C, and the 5700 A line was 
obscured by a cosmic ray in Pal 1-C. The three Pal 1 
[Cu/Fe] ratios agree with the field stars. However, the 
[Cu/Fe] values are not within IS of Ter 7 and Pal 12. 
M67-141's [Cu/Fe] is a bit high, but its 16 errors place it 
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Fig. 13. — [X/Fe] ratios of the elements Cu and Zn versus metal- 
licity. The points are as in Figure [8] 



with the Galactic field stars. 

Three intermediate strength Zn I lines lie in our spec- 
tral region, at 4722, 4811, and 6362 A, and yield consis- 
tent abundances when all three are measured in one star. 
None are detectable in Pal l-I or -II, though all three are 
measurable in M67-141, Pal l-III, and Pal 1-C. Similar 
to copper, both M67-141 and Pal 1 have [Zn/Fe] ratios 
that agree with field stars and GCs but are higher than 
Ter 7 and Pal 12. 

6.6. r- and s-Process Elements 

The heavier elements are formed through neutron cap- 
tures onto iron-peak atoms, either via the rapid (r-) pro- 
cess in Type II supernovae, or via the slow (s-) process, 
e.g. in low mass AGB stars. In the Sun, the percent- 
ages of Y, Ba, La, Nd, and Eu that come from the s- 
Tocess are 75%, 8 5%, 72%, 47%, and 3%, respectively 
Burris et al.ll2000|) : thus Eu is an important r-process 
indicator. The [X/Fe] ratios for the elements that are 
primarily due to the s-process (Y, Ba, and La) are shown 
in Figure 1141 while those with larger r-process contribu- 
tions (Nd and Eu) are shown in Figure ITSl 

6.6.1. s-Process Elements 

Five Y II lines are detectable in the Pal 1 spectra, all 
in the blue: 4884, 4900, 5087, 5200, and 5403 A. The 
first four lines are of intermediate-strength in the Pal 1 
stars while the last is weak. All five lines yield consis- 
tent abundances when measured in one star. Spectrum 
syntheses were performed on the 5087 and 5200 A lines 
for Pal l-I and on the 5200 and 5403 A lines for -II; 
only upper limits are given. While M67-141 agrees well 
with the field stars, the Pal 1 stars are deficient in [Y/Fe] 
compared to field stars and globular clusters. The Pal 1 
[Y/Fe] ratios are in better agreement with those from 
Ter 7 and Pal 12. 

Ba II and La II have three and four lines, respectively, 
in our line list. As mentioned in Section |5.2[ spectrum 
syntheses were performed on all the 6142 A lines to in- 
clude the effect of the blended Ni I line. Another Ba II 
line, 4934 A, was quite strong and was not always below 
the 200 mA limit. The remaining Ba II line, 5853 A, is 
of intermediate to strong strength. The four La II lines 
(5302, 5304, 6320, and 6774 A) are all weak in the Pal 
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Fig. 14. — [X/Fe] ratios of the elements Y, Ba, and La versus 
metallicity. The points are as in F igure [S] writh additional Galactic 
abundances from ISimmerer et al.l 1120045 ^ . 



1 stars. Spectrum syntheses were performed on the 6774 
A line to obtain upper limits for Pal l-I and -II. [Ba/Fe] 
and [La/Fe] show a similar trend in Figure [141 both are 
higher in Pal 1 than in the field stars, but both are in 
agreement with 47 Tuc, Ter 7, Pal 12, and the old GCs, 
particularly Be 29. 

High s-proc ess yields have b een observed in Galactic 
GG (e.g. M4. lYong eralll2008[ ). and have typically been 
explained by primordial variations, i.e. the cluster hap- 
pened to be born in a region where all s-process yields 
were particularly high. However, in Pal 1 not all the 
s-process elements are high: Y, a first-peak s-process el- 
ement, shows the opposite trend from Ba and La, which 
are second-peak elements. To compare the relative con- 
tributions of second-peak to first-peak, [Ba/Y] vs. [Fe/H] 
is shown in the top panel of Fig. 1161 Though only lower 
limits are available for Pal l-I and -II, it is still evident 
that Pal 1 has higher [Ba/Y] ratios than the field stars 
and GGs, though it is in good agreement with the ratios 
for Ter 7 and Pal 12. 

6.6.2. r- Process Elements 

Five Nd II lines are observable in M67-141. Of these 
five, only two (5250 and 5320 A) are strong enough to be 
seen in the Pal 1 stars. Only the latter line was detected 
in Pal l-I; no Nd II lines were found in Pal l-II. The 
6645 A Eu II line required spectrum syntheses for all 
stars; again, only upper limits were obtained for Pal l-I 
and -II. 

Figure [15] shows that M67-141 follows the field star 
trend for both Nd and Eu, while Pal 1, though slightly 
high, is still in agreement with the Galactic field stars 
within its 15 errors. Pal I's [Nd/Fe] ratio may be a bit 
low compared to Ter 7 — no comparisons are available for 
the other Galactic GGs. The [Eu/Fe] of the Pal 1 stars is 
in good agreement with Pal 12, but those are both much 
higher than in 47 Tuc and other Galactic GGs and OGs. 

Sin ce 97% of Eu is produced via the r-process in the 
Sun ([Burris et al . 2000). it is primarily an r-process in- 
dicator. Also, since the site of the r-process is believed 
to be Type II supernovae, which also distribute the a- 
elements, then the [Eu/a] ratio should be correlated. 
The middle panel in Figure [T6l shows [Eu/a] vs. [Fe/H]. 
Pal 1 has higher [Eu/a] ratios than the Galactic field 
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Fig. 15. — [X/Fe] ratios of the el eme nts Nd and Eu versus metal- 
licity. The points are as in Figure [Ml 
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Fig. 16.— [Ba/Y], [Ba/Eu], and [En/a] versus metallicity. The 
points arc as in Figure [141 

stars, GCs, and OCs, but is in good agreement with the 
[Eu/a] ratios in Ter 7 and Pal 12. The excess Eu rela- 
tive to the Galactic stars suggests different yields for Eu 
and a elements, and that whatever the source of these 
different yields, it was similar for Pal 1, Ter 7, and Pal 12. 

6.6.3. s-Process vs. r-Process 

The ratio [Ba/Eu] provides a clue of the ratio of the 
s- to r-process contributions. The bottom panel in Fig- 
ure [16] shows that Pal 1 is in good agreement with the 
Galactic field star distribution, suggesting that in Pal 1 
the s-process contributes to the chemical evolution in a 
similar way as in the Galaxy. It is interesting that Pal I's 
[Ba/Eu] ratios are similar to Pal 12 and Be 20, but that 
those are all slightly lower than the rest of the Galactic 
GCs and OCs. 

7. DISCUSSION 

With the derived abundances of twenty one different 
elements, the guiding questions of this analysis can be 
addressed. Ultimately we wish to understand what type 
of cluster Pal 1 is and whether or not it has been accreted 
from a satellite galaxy. The answers to these questions 
may have implications for galaxy and cluster formation. 



7.1. What Kind of a Cluster is Pal 1? 

As discussed in Section [1] Pal 1 has traditionally 
been classified as a GC, primarily because of its loca- 
tion above the Galactic plane and its high concentration 
parameter. However, Pal I's high [Fe/H] and red hor- 
izontal branch make it similar to the bulge/disk GCs 
(jMackev fc Van Den Be rgh 2005), while its age sets it 
apart from all GCs. In fact, its age and [Fe/H] place it 
between the classical definitions of GCs and OCs. These 
unusual characteristics ha ve caused Pal 1 to be labeled 
as a "transitional cluster" ([Ortolani et al.|[T995[ ). 

In addition. Pal I's brightness and size make it an ul- 
tra faint cluster. A graph of absolute magnitude versus 
hal f-light radius (see Figure 1 in Nicdcrstc-Ostholt et al] 
|2010( ) shows that Pal 1 is fainter than typical globular 
clusters and occupies a position that could be an extrap- 
olation of the dwarf galaxy and ultra faint dwarf galaxy 
trend. It is therefore worth investigating Pal I's similar- 
ities to these different types of objects. 

7.1.1. Basic Properties of Pal 1 

We begin by summarizing the basic properties of Pal 
1 that have been determined in this paper. We have 
derived a average metallicity of [Fe/H] = — 0.60±0.01 for 
Pal 1, a r esult that is in agreeme nt with both the CaT 
estimates (iRosenberg et al.lll998bt ) and the isochrone fits 
(|Saraiedini et al.ll2007[ ). We further find that the cluster 
is not a-enhanced. Defining a to be an average of Mg, 
Ca, and Ti, the average [a/Fe] « 0.0. 

A revised age can be derived using our updated, de- 
tailed abundances. Isochr ones from the Dart mouth Stel- 
lar Evolution Database (jDotter et al.l l2008fl have been 
adopted to derive a new age for Pal 1 using the de- 
rived [Fe/H] « —0.6 and [a/Fe] ~ 0.0 (previously val- 
ues of [Fe/H] = -0 .7 and [a/Fe] = 0.2 were used, 
iSarajedini et al.ll2007| ). The same distance modulus and 
reddening are used as for the photometric atmospheric 
parameters: (to — M)o = 15.65 and E{B—V) = 0.20 (see 
Sec tion | 4.2 | for deta i ls); th e photometry is taken from 
the lSaraiediniet al.l (f200l HST Globular Cluster Trea- 
survF^ 

Fits of the iDotter et al.l (|2008) isochrones are shown 
in Figure [U for ages of 4, 5, and 6 Gyr; the best age 
remains at 5 ifc 1 Gyr. The difference s between these fits 
and those of ISarajedini et afl (|2007| ) are very slight, as 
expected since the adopted [Fe/H] and [a/Fe] combina- 
tions correspond to roughly the same [M/H]. Thus, our 
analysis confirms the previous findings that Pal 1 is in- 
deed a young cluster. We further note that none of our 
target red giant stars lie on the isochrone RGBs; this may 
suggest that our targets are actually horizontal branch 
starso 

No conclusive signs of any abundance spreads are de- 
tected among the four Pal 1 stars, suggesting that Pal 1 
is indeed a simple stellar population. Table [TOl shows the 
mean [X/Fe] ratios in Pal 1 and the abundance spreads 
within the cl uster, according to the method outlined in 
iCohenI ()2004| ). Here a is the dispersion about the mean 

^^ pTttp: //www. astro. ufl.edu/$\sim$ata/public_hstgc/databases. html I 
^■^ If the stars are horizontal branch stars, we would expect to 

obtain slightly different surface gravities, due to, e.g., mass loss. 

However, such differences should have a minor (< 0.1 dex) effect 

on logg. 
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abundance in the cluster and cr (obs) is the mean uncer- 
tainty of an element in a single star. The spread ratio is 
then a comparison of the width of the mean distribution 
to the width of a single star's error estimate: large values 
of the spread ratio imply that the star-to-star variations 
are larger than the average uncertainty for an individual 
star, suggesting that the cluster has a genuine abundance 
spread. The Pal 1 stars do not show a significant spread 
for most elements (a spread ra tio > 1.0 would be con- 
sidered significant: lCohenll2004l ). Cr I has a spread ratio 
> 1.0, but abundances are only available for 2 stars. The 
low spread ratio for Na suggests that a Na/0 anticorre- 
lation is not present, although O is available for only one 
star. Thus, Pal 1 does not appear to show any signs of 
star-to-star variations. This is further confirmed by its 
CMD: multiple populations are not evident and a single 

isochrone fits it well, implying that the cluster is coeval. 
Thus, our analysis finds Pal 1 to be a metal-rich, non- 

a-enhanced, young, chemically homogeneous stellar pop- 
ulation. 

7.1.2. Globular or Open Cluster? 

Given that the cluster types are not well-defined it is 
difficult to distinguish between the two types of clusters, 
particularly when examining chemical abundances. How- 
ever, multiple populations have been observed in nearly 
all Galactic GCs, but not in open clusters. Thus, signs 
of the Na/0 and/or the Mg/Al anticorrelations would be 
positive indicators for a GC. 

No signs of the Mg/Al anticorrelation are detected in 
Pal 1, and O abundances are only available for one of the 
four stars. Pal 1-C. As discussed in Section l673l while Na 
is slightly high and O is slightly low in Pal 1-C, there is 
no significant range in the Na abundances, and the four 
stars do not show evidence for an anticor relation. Cer- 
tain models (e.g. iConrov fc SpergeJI 120111 ) suggest that 
the formation of a second population in a GC depends 
on the cluster mass, since more massive clusters can re- 
tain more gas to form a second populatio n — Pal I's cur- 
rent mass, log(M/MQ) « 3.2 ( Nicdcrst e~Ostholt et all 
|2010( ). is below the critical mass for retaining the gas 
to form a second population. It has been known for 
some time that Pal 1 i s in the process of evaporating 
(jRosenberg et al.lll998al) : lNiederste-Ostholt et all (|2010h 
further show that the mass in Pal I's tidal tails is roughly 
equal to its current mass, suggesting that the cluster 
might have been at least twice as massive in the past. 
In the IConrov fc Spergel (|2011| ) context, however. Pal 
I's lack of a second generation implies that it was never 
more massive than log(M/M0) = 4.0, regardless of its 
current size. This restriction is not necessarily unusual 
for a GC: the GC Pal 12 do es not show signs of abun- 
dance spreads (jCohenI l2004f ) and five intermediate age 
(< 2 Gyr) LMC GCs show no evidence for multiple pop- 
ulations in their CMDs. 

Pal I's detailed abundances can also be compared to 
the well-established GCs in the Galaxy. Its metallicity 
alone puts it in an unusual regime: while there are many 
metal-rich GCs, those are usually centrally concentrated 
in the bulge or disk, while Pal 1 is located much further 
from the center and far above the disk {Rqc = 17.2 kpc 
and Z = 3.6 kpc; lHarrislll996L 2010 edition) Most of 

^^ Recall that previous classifications have labeled Pal 1 as a 



the iron-group (Figures [TT] and [T^ , a- (Figures [5] and 
[9]) , and neutron capture elements (Figures [Ml [151 and 
\W\\ are in good agreement between Pal 1 and the Sgr 
clusters; however, the Na, Al, Cu, and Zn abundances 
(Figures [10] and [T3l) are distinctly different. Thus, the 
chemical pattern of Pal 1 is not clearly similar to either 
the metal-rich Galactic GCs or the Sgr GCs. Of the GCs, 
Pal 1 is most similar to Pal 12 and Ter 7, both known to 
be associated with the Sgr dSph. 

There are several old, metal-poor GCs with detailed 
abundances in the literature, though most are more 
metal-rich than Pal 1. Several of these clusters are lo- 
cated in the outer disk, and may have been accreted from 
a dwarf galaxy or formed during a major merging event 
(|Yong et al.ll2005[ ). Be 20 and 29, the two clusters in our 
sample that lie closest to Pal I's metallicity, have similar 
abundance ratios as the Galactic stars, with the excep- 
tions of the s-process neutron capture elements. Thus, 
the only elements that are clearly discrepant between Be 
20/29 and Pal 1 are the a-elements (Figures [8] and [9]) 
and Eu (and by extension [Eu/a] and [Ba/Eu]; Figures 
[T5l and [TBI) ■ C'u, Zn, and Y abundances are not available 
for Be 20 and 29; however, Pal I's Na and Al abundances 
(Figure [To]) agree better with Be 20 and 29 than with the 
Sgr clusters. 

It is then natural to ask if Pal 1 looks more like the GCs 
Pal 12 and Ter 7 or the GCs Be 20 and 29. In terms of 
the a-elements and neutron capture elements Pal 1 looks 
more like the GCs, while Na and Al agree more with the 
GCs. However, chemically comparing Pal 1 to the two 
cluster types is not an accurate way of determining Pal 
I's cluster type, as Pal 1, Be 20, and Be 29 are not con- 
clusively associated with dwarf galaxies while Pal 12 and 
Ter 7 are associated with the Sgr dSph. Furthermore, 
the chemical abundances of stars depend on their host 
galaxies' star formation histories, etc., and these quan- 
tities are unique to individual galaxies. Thus, unless all 
of the clusters originated in the Sgr dSph, this compari- 
son is not straightforward. We leave the question of Pal 
I's origin to Section [721 and we therefore focus on other 
parameters that may distinguish GCs from OCs. 

Table [11] shows various quantities that can be com- 
pared between the clusters: the Galactocentric radius, 
Rqc (in kpc); the distance above the Galactic plane, Z 
(in kpc); distance from the Sun, d (in kpc); the age (in 
Gyr; though note that these ages are assembled from 
different sources, and may not be valid for comparisons); 
metallicity, [Fe/H]; absolute visual magnitude. My; half- 
light radii, rt (in pc; note that the distance was used to 
convert the half-light radius from arcminutes to pc); and 
concentration parameter, c. The last two parameters are 
only available for the GCs. Note that w e also include 
Whiting 1, a GC associated with Sgr (jCarraro et al.l 
(2007), though there are no chemical abundances avail- 
able for this cluster at this time. The values of Rqc are 
similar for all the clusters in the table. The OCs, while 
more than 1 kpc away from the plane, are still closer 
than the GCs, with Pal 1 lying between the two. The 
ages and metallicities are similar, but again Pal 1 is still 
quite young for a GC. Though Pal 1 is fainter than the 



bulge/disk cluster (e.g. IMackev & Van Den Berghll2005l ') based on 
its [Fe/H] and HB morphology rather than its location, but its 
unusual location does distinguish it from the bulge/disk GCs. 
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TABLE 10 

Mean abundances and spread ratios for th e four stars in Pal 1, 

ACCORDING TO THE METHOD Of ICOHEnI II2004I) . A SPREAD RATIO > 1 

INDICATES THAT A SIGNIFICANT SPREAD EXISTS; PAL 1 SHOWS NO 

EVIDENCE FOR ANY ABUNDANCE SPREADS. SEE THE TEXT FOR A 

DESCRIPTION OF HOW THESE VALUES ARE CALCULATED. 



Mean [X/Fc]^ Number of stars 



(T(obs) Spread Ratio 



Fe I 


-0.60 


4 


0.01 


0.06 


0.26 


Fell 


-0.60 


4 


0.10 


0.21 


0.46 


Na 


0.21 


4 


0.04 


0.21 


0.21 


Mg 


-0.06 


4 


0.07 


0.16 


0.44 


Al 


0.04 


2 


0.05 


0.10 


0.49 


Si 


0.16 


4 


0.07 


0.15 


0.45 


Ca 


0.10 


4 


0.09 


0.13 


0.75 


Sc 


0.25 


4 


0.05 


0.19 


0.25 


Til 


-0.06 


4 


0.07 


0.21 


0.33 


Tin 


-0.04 


3 


0.11 


0.21 


0.52 


V 


0.06 


4 


0.08 


0.12 


0.68 


Crl 


-0.09 


2 


0.12 


0.10 


1.20 


Cm 


0.18 


2 


0.01 


0.21 


0.07 


Mn 


-0.16 


3 


0.05 


0.18 


0.27 


Co 


-0.04 


2 


0.09 


0.20 


0.47 


Ni 


0.06 


4 


0.03 


0.12 


0.22 


Cu 


-0.03 


3 


0.06 


0.28 


0.23 


Zn 


0.12 


2 


0.05 


0.15 


0.34 


Y 


-0.41 


2 


0.06 


0.12 


0.53 


Ba 


0.24 


4 


0.04 


0.20 


0.18 


La 


0.33 


2 


0.13 


0.15 


0.85 


Nd 


0.15 


3 


0.04 


0.32 


0.12 


Eu 


0.50 


2 


0.0 


0.20 


0.0 



^ [X/H] is given instead for Fc I and Fc ll. 

OC Be 29, its brightness is comparable to the GC Whit- 
ing 1. Finally, Pal 1 has a very small half-light radius 
compared to Pal 12 and Ter 7, but in agreement with 
Whiting 1, and Pal 1 is very concentrated, much more 
SO than Ter 7 and Whiting 1. 



The distinctions between OC and GC in Table [Tl] 
seem to be based primarily on distance from the Galac- 
tic plane. Considering that Pal 12, Ter 7, and Whiting 
1 have been accreted by Sgr, this is not a particularly 
valid criterion. Consequently, the distance above the Sgr 
plane, |^Sgr|, is also plotted for clusters that are (or might 
be) a ssociated with the Sgr dSph (from lLaw fc Maiewskil 
I2010t note that Pal I's |^sgr| is also included, although 
it is most likely not associated with Sgr) . Examined this 
way, the clusters look very similar, i.e. the young, metal- 
rich GCs that are close to the Sgr plane are similar to the 
old, metal-poor OCs that are close to the Galactic plane, 
suggesting that these clusters may he classified as open 
(or inte rmediate-aged) cluste rs in the Sgr frame of ref- 
erence. iCarretta et al.l ()201C1[ ) also cast doubt upon the 
GC classification for Pal 12 and Ter 7, given the lack of 
a definite Na/0 anticorrelation. Under the Carretta et 
al. definition we are forced to conclude that Pal 1 is not 
a bona fide GC; however, it is not an obvious Galactic 
OC either, and therefore remains as an unusual cluster. 

7.1.3. A Connection with Ultra-Faint Dwarfs? 

Given Pal I's identification as an ultra faint cluster 
and its proximity to the ultra-faint dwarfs (UFDs) in a 
plot of a bsolute magnitude vs. half- light radius (Fig- 
ure 1 in .Niederste-Ostholt et al.ll2010D it is tempting to 
consider whether Pal 1 is related to these other ultra- 
faint objects. The UFDs are distinguished by their 
larger-than-expected velocity dispersions, which imply 



high mass/light ratios and large amounts of dark mat- 
ter. The velocities for the five Pal 1 stars (including 
Pal 1-IV) listed in Table [5] imply a velocity dispersion of 
3.6 lb 1.5 km s~^ (using the formula from IWalker et aTl 
I2006f ). assuming that all five velocities are orbital veloci- 
ties of cluster members. Under this assumption. Pal I's 
velocity dispersion is in agreement with the velocity dis- 
persions of several UFDs, including Segue II, Leo V, Leo 
IV , and Hercules (see the sum mary of UFD properties 
bv lMcConnachie fc Cot^lMoh . 

As tempting as this comparison may be, it is not al- 
together appropriate. Firstly, we now know that Pal 1 
has tida l tails and is therefore not in dynamical equi- 
librium (jNiederste-Ostholt et alj [20100 . Secondly, only 
one star (Pal 1-IV) has a discrepant velocity, and this 
difference is easily explaine d if Pal 1-IV is a binary; 
iMcC onnachie fc Cotel (j2010f) showed that the presence 
of binaries can significantly boost the velocity disper- 
sion of a low-mass system. Alternatively, Pal 1-IV could 
be a nonmember — its location in the CMD (Figure [T]) 
is separate from the other stars — however its consistent 
chemistry and very close radial velocity make this seem 
unlikely. Therefore, Pal 1 remains an ultra faint cluster, 
and though it would be interesting if it were an extrapo- 
lation of the UFDs, its current velocity dispersion is best 
explained by its tidal disruption or by the presence of a 
binary. 

7.2. Is Pal 1 Extragalactic? 

Pal 1 has been tentatively associated wi th both the 
Galactic Anticenter Stellar Stream (GASS) (jCrane et al.l 
200_3') and the C anis Major (CMa) ove rdensity (e.g. 
iMartin et al.ll2004 iForbes fc BridgeslboiOD . two streams 
of stars that could be from disrupted dwarf galaxies. If 
it is associated with either of these systems, then Pal 1 
should have similar chemical abundances to their associ- 
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TABLE 11 
Various parameters for the GCs and OCs used in this s tudy. If no re ferences are given then the 

DATA comes from THIS STUDY (iN THE CASE OF PAL 1) OR THE [HARRIS (1993) CATALOG (2010 EDITION). NOTE 

THAT THE AGES ARE ASSEMBLED FROM DIFFERENT SOURCES, AND THER EFORE MIGHT NO T BE DIRECTLY 

COMPARABLE. AlL VALUES OF {Zsgrl ARE FROM ILAW fc MAJEWSKJ I|2010I) . 



Cluster Traditional Rqc Z d age [Fe/H] My r^"" 

Classification (kpc) (kpc) (kpc) (Gyr) (pc) 



I^Sgrl References'^ 
(kpc) 



Pal 1 

Pal 12 
Ter 7 
Whiting 1 
Be 20 
Be 29 
Saurer 1 



GC 
GC 
GC 
OC 
OC 
OC 



17.2 
15.8 
15.6 
34.5 
15.8 
21.6 
19.2 



3.6 
-14.1 

-7.8 
-26.3 
-2.5 
1.8 
1.7 



14.2 
19.0 
22.8 
30.1 



5.0 
9.5 
8.0 
6.5 

5.8 
4.5 
5.0 



-0.6 
-0.80 
-0.60 
-0.65 
-0.44 
-0.54 
-0.38 



-2.52 
-4.47 
-5.01 
-2.46 
-2.06 
-4.64 



1.8 
9.1 
4.9 
1.8 



2.57 
2.99 
0.93 
0.55 



9.8 
3.24 
1.09 
0.215 

1.79 
5.25 



1, 2 
1, 3 

4 

5, 6, 7 

4, 6, 8 

4, 9 



^ The distance d was used to convert the half-light radius (rh) from arcminutes to parccs. 

^ References. (l)lDottcr ct al. (2010): (2) Cohcnl ll2004f) : (3) Sbordonc ct al. (200 5^; (4) from th e sources as sembled by 
[Forbes & Bridges (2010); (5) Andrcuzz i ct al.. L20lB : (6) lYong ct al, ( 2005. ): (7) .De Marchi et al.l 1 D0O&') : (8) .Lata et al.l 



1I2OO2I) : (9) IFrinchabov et al.l 120061) 



ated field stars and clusters and to the stars in other low 
mass systems. 

Figures [T7] through [^S] compare the Pal 1 abundances 
to those of stars in the nearby dwarf galaxies. Note 
that only the average Pal 1 abundance (red star symbol) 
is plotted, with the mean abundance error for [X/Fe] 
(e.g. (T(obs) in Table [T0|) . The Galactic field stars 
are shown in grey. Dwarf galaxies are also included, 
with their field stars as small symbols and their GCs 
as larger open symbols. Again, abundances have been 
averaged over an entire clust er. Black points are For- 
nax stars (jLetarte et al.|[20Tol ): cyan triangles are Sagit- 
tariu s (Sgr) stars (field s tars are from [S bordon e et al.l 
120071 iMonaco et aLll2007l andlChou et al . 2010a. while 
GC st ars are from iCohenI 120041 and ISbordon e et al.l 
|2005b); blue squares are Large Magellanic Cloud stars 
(fie ld stars from Pompcia et al. 2008; cluster stars from 
[Johnson et al.l l2006l and iMucciarelli et al.l I2008D ; green 
crosses are Galact ic Anticenter Stellar Stream (CASS; 
iChou et "ani2010bl ) stars; and orange di amonds are field 
stars in the Canis Major overdensity (jSbordone et al.l 
l2005al) . While other dwarf galaxy data are available in 
this metallicity range, plotting all of it makes these plots 
nearly illegible, and we have therefore included only the 
most relevant dwarf data (i.e. galaxies with many stars 
and many elemental abundances in the same metallicity 
range as Pal 1) for this discussion. 

7.2.1. General Comparisons: Dwarf Galaxy or MW cluster? 

In terms of the a-elements (Figures [T71 and [T5)) . Pal I's 
slightly low [a/Fe] ratios place it with the dwarf galax- 
ies (with the possible exception of Ca). This suggests 
that Type la supernovae began to contribute to Pal I's 
host environment at a lower [Fe/H] than for stars in the 
Galaxy. In general, Pal I's [a/Fe] ratios agree with the 
Sgr and LMC clusters. The slightly subsolar [Ti/Fe] 
value agrees with the CASS stars while the slightly low 
[Mg/Fe] and [Ti/Fe] and slightly higher [Ca/Fe] and 
[Si/Fe] ratios agree with the [Fe/H] = -0.44 CMa star. 
The latter CMa star is also similar to Pal 1 in Na and 
Al (Figure [TO]) : for the few available [Al/Fe] abundances 
Pal 1 is higher than the Sgr and LMC clusters. 

Of the iron-peak elements (Figures |20| and |2T|) . Pal I's 
[Sc/Fe] remains clearly distinct from the dwarf galaxy 
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Fig. 17. — [X/Fe] ratios of «-elements O, Mg, and Si versus 
metallicity for the stars in our study as well as Galactic field 
stars (grey; as in Figure [Sjl an d stars in ot her galaxies. Black 
points are Fornax stars jLetarte et al.l [20101 ): cyan triangles are 
Sagittarius (Sgr) st ars (field stars are fro m "Sbordonc ct al. 200^, 



Monaco et al.r2007|, andiChou et al.l2010al . while GC stars arc from 
Cohenll2004l and ISbordone et al. 2005b): blue squares arc Large 
Magellanic C loud stars (field stars from Pompcia ct al. 2008; clus- 
ter stars from "Johnson et al.ll2006l and lMucciarelli et al.ll2008. ); and 
orange diamonds are field stars in the Canis Major overdensity 
|[?bordone et al.ii2005al ). 

stars. Pal I's [V/Fe] ratio is in good agreement with 
the dwarfs, with the exception of the CMa stars, which 
may be due to HFS corrections. Overall the Pal 1 Fe- 
peak elements show the best agreement with the Galac- 
tic stars rather than the dwarf galaxy stars. Variations 
in [Fe-peak/Fe] with metallicity have been suggested to 
be due to metallicity-dependent supernovae yields, e.g. 
metal-poor Typ e la supernovae will produce less Mn than 
metal-rich ones (|Cescutti et al.ll2008:) . as seen in the Sgr 
field stars. However, Pal I's agreement with the Galactic 
and dwarf galaxy stars at a similar metallicity suggests 
no significant dispersion in the Type la contributions. 
The more peculiar elements, Cu and Zn (Figure [55]), are 
in better agreement with the Galactic stars and possi- 
bly CMa, and are clearly distinct from the Sgr and LMC 
stars. 

Finally, Pal I's neutron capture elements (Figures [231 
l24l and |25|) agree better with the dwarf galaxies than 
with the Galactic stars. Pal I's low [Y/Fe] and high 
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Fig. 18. — [X/Fe] ratios of the Q-eleinents Ca and Ti, as well as 
of the average a from Mg, Ca, and Ti. The Ti abundance is an 
average of Ti I and Ti II, except for Pal l-I, whose Ti I abundance 
is shown since i ts T i II abundance is abnormally high. The points 
are as in Figure lTTl with the addition of G alactic Anticenter S tellar 
Stream (GASS; green crosses) stars from lChou et al] l|2010bl V 
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Fig. 19. — [X/Fe] ratios of the light elements Na and Al versus 
metallicity. The points are as in Figure [iTl 
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Fig. 20. — [X/Fe] ratios of the iron-peak ele ment s Sc, V, and Mn 
versus metallicity. The points are as in Figure lTTl 

[Ba/Fe], [La/Fe], and [Ba/Y] ratios suggest that like the 
Sgr dSph, Pal 1 was enriched by metal-poor AGB stars. 
(This effect is also seen in the [La/Y] ra tios of the Sgr 
stars; iMcWilliam fc Smecker-Hand 120051) . With an ex- 
cess of neutrons per Fe atom in metal-poor AGB stars, 
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Fig. 21. — [X/Fe] ratios of the iron-peak elements Co, Cr, and 
Ni versus metallicity. Both Cr I and Cr II a re s hown (Cr I in red, 
Cr II in yellow). The points are as in Figure [Tt] 
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Fig. 22. — [X/Fe] ratios of the el eme nts Cu and Zn versus metal- 
licity. The points are as in Figure [TtI 
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Fig. 23. — [X/Fe] ratios of the elem ents Y, Ba, and La versus 
metallicity. The points are as in Figure [TSl 



more heavy second-peak s-process elements such as Ba 
and La can be cre ated, leaving a deficit of first-peak e l- 
ements like Y (Gal lino et al.lll998l : iBisterzo et al.ll2010( ). 
This further suggests that Pal 1 has not been enriched 
by its own metal-rich AGB stars, which is perhaps not 
surprising given its lack of evolved stars (see Figure 1). 
The GASS field stars do not have high [La/Fe] hke Pal 1, 
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Fig. 24. — [X/Fe] ratios of the el eme nts Nd and Eu versus metal- 
licity. The points are as in Figure (TT] 
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Fig. 25. — [Ba/Y], [Ba/Eu], and [Eu/q] versus metaUicity. The 
dashed red Une shows the r-process-only [B a/E u] ratio (from 
IBurris et al.ll2000l ). The points are as in Figure [iT] 

suggesting that it does not have an excess of second-peak 
to first-peak s-process elements. The CMa stars do seem 
to have high [Ba/Fe] and [La/Fe], and arguably have low 
[Y/Fe], in agreement with Pal 1. 

Pal 1 also has shghtly higher [Eu/Fe], [Nd/Fe], and 
[Eu/of] ratios than Galactic stars, in agreement with the 
LMC, Sgr, and Fornax. While the CMa stars do show an 
excess of Nd, they do not appear to have an excess of Eu, 
and CMa therefore has a normal value for [Eu/a]. This 
suggests that dwarf galaxies have an additional source 
of r-process elemen ts compared to the Galaxy (see e.g. 
iLetarte et al.|[20Tol ). This is further supported by Pal 
1 and the dwarf galaxies' high [Eu/a] ratio; the stan- 
dard model suggests that massive Type II supernovae 
are responsible for dispersing the a-elements, while lower 
mass (8-10 Mq) Type II supernovae are the sole site of 
r-process elements like Eu. If the dwarf galaxies had 
similar conditions to the Galaxy, then the [Eu/a] ratio 
should remain flat. 

Despite its high values for [Ba/Fe] and [Eu/Fe], the 
[Ba/Eu] ratio in Pal 1 agrees better with the Galactic 
stars than with the dwarf galaxies, though the Fornax 
dispersion and Pal 12 do extend down to Pal I's value. 
Again CMa's [Ba/Eu] values are in fair agreement with 
Pal 1. 

As discussed earlier, at solar metaUicity Ba is primarily 



an s-proc ess element , while Eu is primarily an r-process 
element (jBurris et al.ll2000[ ). The ratio [Ba/Eu] should 
therefore provide an indication of the relative contri- 
butions of the s-process to the r-process. Before the 
s-process begins to contribute Ba is initially produced 
solely through the r-process; as time goes on the s-process 
creates more Ba but little Eu, and [Ba/Eu] increases. 
The theoretical r-process-only lower limit for [Ba/Eu] is 
shown in Figure[25l Pal I's slightly high Ba and Eu abun- 
dances lead to a [Ba/Eu] value that agrees with the clas- 
sical prediction for a system that has had little contribu- 
tion from the s-process. Clearly this cannot be the case: 
Pal I's [a/Fe] ratios show that Type la supernovae have 
contributed to the system, and therefore AGB stars must 
have also contributed since the timescale for evolution of 
Type la supernovae and AGB stars is th ought to be sim- 
ilar in chemical evolution models (e.g. [T ravaglio et al.l 
[200l IMatteucci et"al] [2001 IZolotov et al. 2010 ). The 
[Ba/Y] ratio further implies that metal-poor AGB stars 
must have contributed to the system. Thus, the high 
[Eu/Fe] and [Ba/Fe] together with the high [Eu/a] sug- 
gest that there is an extra source of r-process elements 
in the Pal 1 system, beyond what is present for typi- 
cal Galactic stars. It is unclear where the additional r- 
process abundances come from; some possibilities include 
shot noise (i.e. simple inhomogeneous mixing of a nearby 
supernova with a progenitor in the 8-10 M© range), an 
unusual star formation history that may have included 
an effectively truncated IMF (missing the most massive 
stars that contribute significantly to the a elements), or 
even variable nucleosynthetic yields. However, each of 
these possibilities has other consequences for the chemi- 
cal abundance ratios that are not clearly indicated in Pal 
1. Whatever the case. Pal 1 has clearly not followed the 
"classical" chemical evolution models. 

Based on these chemical comparisons, we conclude 
that Pal 1 likely originated in a dwarf satellite that 
was later accreted by the Milky Way (about 500 Myr 
ago: iNiederste-Ostho'lt et al.ll2010D . Though the light el- 
ements {Z < 30) are arguably in agreement with Galac- 
tic disk stars, the neutron capture elements are distinct 
enough for us to conclude that Pal 1 does have an extra- 
galactic origin. 



7.2.2. Comparisons with Specific Galaxies: 

Galaxy? 



Pal 1 's Host 



We now attempt to chemically link Pal 1 to its host 
galaxy, as has been done with e.g . Pal 12 (jCohenll2004D 
and Ter 7 ([Sbordone et al.|[2005bD . In particular we ex- 
amine the chemical abundances of the CMa overdensity, 
the GASS, and the LMC intermediate-age clusters. 

Given the low [La/Fe] abundance in the GASS, it is 
not likely that Pal 1 is a member of this stream. The few 
stars analyzed in the CMa overdensity are more promis- 
ing (even if the stars analyzed so far in CMa are more 
metal- rich): the general trends in the a-abundances, the 
odd-Z elements Na, Al, Cu, and Zn, the iron-peak el- 
ements (with the exception of Sc and V possibly due 
to HFS corrections), and the neutron capture elements 
(possibly with the exception of Eu) are in fair agreement. 
Based on this analysis it is possible that Pa l 1 could be as- 
sociat ed wi th this stream, as suggest ed by iSaviane et al.l 
(|20T0l) and iForbes fc Bridges! (|2010l) . The slight differ- 
ences between individual elements could be due to the 
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choices of spectral lines and atomic data. However, it 
must be noted that none of these stars are guaranteed 
memb ers of CMa, as disk co ntamination is extremely 
likely (jSbordone et al.l l2005al ). In particular, the star 
with an [Fe/H] closest to Pal 1 seems identical to the 
Galactic disk stars, and may not be a true member of 
the stream. Analyses of more stars in the CMa overden- 
sity, particularly in the metallicity range of Pal 1, would 
help establish potential membership. 

A comparison with the LMC intermediate-age clus- 
ters must also be considered, as they are also young 
and metal-poor. In addition, there are several low- 
mass GCs that appear to be similar to P al 1: in their 
sampl e of sixteen intermediate-age clusters. iMilone et al.l 
(|2009f) found five with no signs of multiple populations 
in their CMDs, suggesting that they were not massive 
enough to retain gas fo r a second generation of stars 
(|Conrov fc Spergellboill see Section [LT|) . Though the 
Na, Al, Cu, and Zn abundances in these clusters are 
slightly different from Pal 1 , the rest of the elements are 
in good agreement. Pal 1 does not appear to have been 
accreted from the LMC; however, it is possible that its 
formation mechanisms were similar to the intermediate- 
age LMC clusters, i.e. that Pal 1 was accreted from a 
dwarf galaxy similar to the LMC. 

Finally, with the exceptions of [Na/Fe], [Al/Fe], 
[Cu/Fe], and [Zn/Fe] (which are all higher in Pal 1), Pal 
1 agrees well with Pal 12 and Ter 7 in Sgr. Thus, regard- 
less of its association with either the CASS or CMa, we 
conclude that Pal 1 likely originated in a fairly massive 
dwarf satellite (i.e. a satellite that had a mass somewhere 
between Sgr and the LMC). 

8. SUMMARY AND CONCLUSIONS 

Detailed chemical abundances have been found for 
twenty one elements in four red giant stars in the un- 
usual GC Palomar 1. Our findings are summarized as 
follows: 

• Pal 1 is an outer halo, young (5 ± 1 Gyr) cluster 
whose age, metallicity, location, and structural pa- 
rameters distinguish it from the standard Galactic 
globular or open clusters. However, Pal 1 does ap- 
pear to be similar to the known extragalactic GCs 
and GCs that have been accreted during mergers 
(i.e. from the Sgr dSph). 

• Pal I's young age and single stellar population sug- 
gest a resemblance to the low-mass intermediate- 
age LMC clusters. Chemical abundances further 
suggest that Pal 1 may have shared a similar for- 
mation history. 

• Pal 1 shows several unusual chemical characteris- 
tics, including 

— Lower [a/Fe] ratios than Galactic stars of the 
same [Fe/H] 



— A lack of evidence for a Na/0 anticorrelation, 
though marginally high Na is found 

— Similar [Fe-peak/Fe] ratios to Galactic stars 

— An excess of second-peak s-process neutron 
capture elements over first-peak elements 

— Low [Ba/Eu] and high [Eu/a] values that sug- 
gest Pal I's host galaxy had an additional r- 
process site (possibilities could be inhomoge- 
neous mixing, an effectively truncated IMF, 
i.e. one that is missing the most massive stars, 
or variable nucleosynthetic yields). 

• Chemically, Pal 1 does not behave like the typ- 
ical Galactic bulge/disk GCs or the old, metal- 
poor GCs. The closest agreement seems to be with 
the Sgr clusters Pal 12 and Ter 7 and the LMC 
intermediate-age clusters; however the Na, Al, Cu, 
and Zn abundances do not agree with the LMC or 
Sgr clusters. 

• Comparing the [X/Fe] of the Pal 1 stars to those 
of stars in known streams show that 

— It is unlikely that Pal 1 originated in the 
Galactic Anticenter Stellar Stream, given the 
differing [La/Fe] ratios 

— Pal 1 may have originated in the Canis Ma- 
jor overdensity, if the [X/Fe] ratios of the 
three possible CMa stars are extrapolated to 
slightly more metal-poor stars. 

Overall, we conclude that Pal 1 likely had an extra- 
galactic origin, though its chemistry remains unique com- 
pared to the known globular and open clusters. 
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